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SUMMARY 


Laser  light  conducted  by  optical  fibers  in  a  liquid  medium  is  modulated  by  acoustic  pressures. 
This  phenomenon  is  exploited  to  make  an  acoustic  sensor  of  underwater  sound  by  demodulating  the 
light  through  a  heterodyning  procedure.  In  application  to  underwater  detection  of  ocean  acoustic  sig¬ 
nals  the  problem  arises  as  to  the  response  of  bare  fibers  to  fluctuations  of  ocean  temperature.  Quantita¬ 
tive  estimates  are  necessary  of  the  extreme  case  where  the  fibers  are  bare,  the  temperature  is  instan¬ 
taneously  sensed  throughout  the  fiber,  and  statistical  equilibrium  is  present.  The  effects  of  conductivity 
are  deliberately  omitted. 

The  procedure  for  making  these  estimates  followed  these  steps:  First,  a  derivation  was  made  of  a 
formula  for  the  signal/noise  ratio  of  a  generic  fiber-optic  hydrophone  which  included  all  noise  contribu¬ 
tions,  shot  noise,  amplifier  noise,  Johnson  noise,  laser  jitter  noise,  environmental  temperature  noise, 
and  environmental  mechanical  stress  noise;  this  report  calculates  numerical  estimates  of  the  last  two 
types  of  noise  to  show  their  phase-distorting  effects.  Second,  the  bare  fiber-optic  sensor,  in  the  form  of 
a  single  length  of  fiber,  was  reviewed  in  terms  of  its  instrument  temporal  and  spatial  windowing. 
Numerical  estimates  are  shown  to  be  strongly  dependent  on  how  much  of  the  wavenumber— frequency 
power  spectrum  of  internal  waves  is  accepted  by  these  windows.  Third,  the  temperature  fields  of  the 
ocean  were  categorized  by  adoption  of  pertinent  Garrett  and  Munk  spectra  of  internal  waves,  including 
the  effects  of  tow  speed.  It  is  shown  from  these  spectra  that  for  ordinary  tow  speeds  the  ocean  tem¬ 
perature  field  is  very  nearly  frozen,  except  for  the  very  shortest  and  very  longest  wavelengths.  Fourth, 
for  specific  lengths  of  bare  fibers  at  specific  depths  quantitative  estimates  were  made  of  the  fiber 
response  to  ocean  temperature  fluctuations.  The  results  are  shown  in  tabular  form.  Briefly,  the  rms 
phase  fluctuation  for  fibers  at  100-m  depth  is  estimated  to  be  4.04  rad  for  a  1000-m  towed  fiber,  3.83 
rad  for  a  100-m  dropped  (i.e.,  vertical)  fiber,  and  1.6  to  19.6  rad  for  a  moored  fiber,  depending  on  the 
duration  of  observation  of  internal  waves.  Fifth,  a  review  was  made  of  microtemperature  scales  in  the 
ocean  and  quantitative  estimates  of  their  effects  on  phase  distortion.  It  is  shown  that  a  fiber  0.4  cm 
long  at  a  depth  of  30  to  60  m  exhibits  a  phase  distortion  of  1.7  rad,  while  a  0.01-cm  fiber  exhibits  a 
0.29-rad  distortion,  all  attributed  to  microtemperature  fluctuations  on  bare  fibers. 
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INTRODUCTION 

A  fiber-optic  hydrophone  based  on  heterodyne  detection  consists  of  three  essential  parts:  a  laser, 
two  fibers,  and  a  photodetector.  At  the  output  of  the  laser  the  beam  is  split  in  two,  and  each  beam  is 
coupled  to  a  fiber:  one  fiber  carries  the  reference  becm\  the  second  fiber  is  exposed  to  the  acoustic 
wavefield,  which  modulates  the  laser  light  to  form  the  signal  beam.  The  two  fibers  feed  simultaneously 
into  the  photodetector,  which  develops  an  electric  current  component  proportional  to  the  phase 
difference  between  the  two  beams.  This  signal  current  is  then  amplified  by  a  preamplifier,  which 
delivers  a  voltage  proportional  to  the  acoustic  signal. 

One  problem  in  the  analysis  of  a  fiber-optic  hydrophone  is  to  determine  its  signal-to- noise  ( S/N ) 
ratio,  which  is  a  measure  of  its  acoustic  performance.  A  second  problem  is  to  obtain  quantitative  esti¬ 
mates  of  the  noise  induced  in  the  signal  beam  by  environmental  disturbances  such  as  temperature  fluctua¬ 
tions  and  mechanical  stresses.  Both  of  these  problem*  are  investigated  in  this  report. 

PROCEDURE  OF  ANALYSIS 

An  optical  fiber  of  length  L  carrying  the  signal  beam  is  assumed  to  be  disturbed  by  environmental 
factors.  It  is  desired  to  calculate  the  effect  of  these  disturbances  on  the  stability  of  the  phase  of  the  sig¬ 
nal  beam  and  from  this  to  determine  the  S/N  ratio  of  the  fiber.  The  procedure  of  analysis  used  in  this 
report  to  make  these  calculations  will  be  based  on  autocorrelation  functions  and  power  spectra  of  the 
random  environmental  disturbances.  To  illustrate  this  procedure,  let  the  disturbance  be  a  random  tem¬ 
perature  field  80 (x,  t).  The  corresponding  random  phase  A0(x,  t )  of  the  signal  caused  by  this  field  is 
derived  from  it  in  the  following  way. 

Relation  Between  Mean-Square  Phase  Fluctuation  and  Mean-Square  Temperature  Fluctuation 

Let  5</>0)  be  the  j  th  realization  of  the  phase  fluctuation  over  a  length  L  of  a  fiber  due  to  tempera¬ 
ture  fluctuation  at  a  point  x  of  the  fiber.  Then 

8^o>  *  r  c«(x)  d*' 

where  C#  is  a  temperature-to-phase  conversion  factor,  with  units  (rad/°C)/m.  Assume  O(x)  **  Cw  = 
const  and  define  a  quantity  X  such  that 

=  C„  X(>\  XU)  =  J0  dx. 

Now  the  mean-squared  value  of8<£^  is 

<b<t>}  >  =  C„2  J*o  <h9(x)b9(x’)>  dxdx'. 

But  by  the  definition  of  autocorrelation 

<b0(x)bO(x’)>  =  <502>  pw(x  -  x'), 
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in  which  <802>  is  the  variance  and  p&y(x  -  jc')  is  the  nondimensional  autocorrelation  of  the  random 
function  5 0(x).  Now  following  Taylor  [1]  one  can  write 

r  <80(x)  89(x')>  dx  =  <80(x')X>  =  \  -fr  <X2>. 

2  dx 

Since  we  shall  always  take  p  to  be  an  even  function,  we  write p(x  —  x')dx  =  p(f  .  Then 

<X2>  =  2<802>J’O  f*p({)d{dx\ 

Two  limiting  forms  are  immediately  available: 

•  If  L  is  short  enough  for  p(f)  to  be  about  unity,  then 

<X2>  =  2<80[>  f0  f*  d(  dx '  =  <89[>L2. 

This  is  the  coherent  case. 


•  If  L  is  so  long  that  many  uncorrelated  inhomogeneities  of  temperature  are  present  over  the 
fiber,  then  there  is  an  integral  scale  Lc  such  that 

<*2>  =  2 <80l>  f0  dx'  —  2 <802l>LcL. 


This  is  the  incoherent  case.  Between  these  extremes  we  can  write  another  form,  obtainable  from  the 
above  forms  by  interchange  of  order  of  integration  [2]: 


<X2>  =  2  <80[>  tj0 


_  jL 

L 


In  the  following  applications  we  shall  be  mainly  interested  in  cases  where  L  »  Lc. 
phase  fluctuation  due  to  temperature  fluctuation  per  unit  of  length  L  then  becomes 


V<S0L> 


=  C„L  -V /2y-  J<802L>,  L, 


<<L 


The  mean-square 


In  the  literature  of  temperature  statistics  in  the  ocean,  the  integral  scale  Lc  of  temperature  fluctuations 
is  only  sparsely  known.  Since  the  formula  for  <8<f>l>  is  valid  for  Lc<<  L,  we  shall  have  to  resort  to 
several  approaches  to  overcome  the  limited  information  about  Lc : 


•  Wherever  it  is  known,  we  use  Lc. 


•  When  Lc  is  not  known,  we  can  assume  a  special  case, 

ILL,  -  0(1  )m2,  Lc  «  L 
V  <8  <f>l>  —  C9LJ<b02L>. 


When  Lc  is  not  known,  we  assume  that  the  autocorrelation  is  Gaussian,  p(f)  =  e 


-t2Hho. 


then 

<8<f>l  >  =  Ch  <802>  [2LLm  -2  Ll  (1  -  e~L*nth) ,/J. 

The  correlation  length  L must  be  assumed.  (See  Appendix  A.) 

We  prefer  the  second  approach  and  shall  use  it  in  calculations  below.  A  similar  derivation  may  be 
made  for  the  case  of  rms  phase  fluctuation  due  to  time-varying  temperature  fluctuations  for  nearly 
incoherent  summation, 
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•\J 8<t> 

T 

11 

“H 

C7fc  / — 

wheie 

—  conversion  constant 

(units: 

Te 

=  correlation  time  of  86 (t) 

(units: 

s), 

T 

=  integration  time 

(units: 

s),  and 

<86>f> 

=  variance  of  86 (/) 

(units: 

°C2). 

When  both  temporal  and  spatial  effects  are  considered  together,  the  combined  variance  <86 >[T>  must 
be  used.  The  result  is 

-  c,tr 

where  CvLr  =  conversion  constant,  in  (rad/°C)/(m/s).  This  formula  applies  only  when  Lc  «  L  and 
Tc  <<  T.  Again,  when  Tc  is  not  known  we  may  choose 

2  7T,  =  0(l)s2,  Tc  «  T, 

so  that 

•\J  <8<f>  f>  =  Cf)  7*\/  <86g>. 

A  more  drastic  assumption  is 

ILL,  TTC  —  0(1)  m2-  s2, 

so  that 

\! <&<t> lr>  ~  CHn  yj <bO[r> . 

It  is  to  be  noted  that  rms  phase  fluctuations  can  be  reported  per  unit  length,  per  unit  time,  or  per  unit 
length  per  unit  time  only  if  the  correlation  length  and  correlation  time  are  known.  If  the  assumptions 
are  made  that  ILL,  =  0(l)m2  and  2TT,  =  0(1  )s2,  then  the  rms  phase  fluctuation  refers  to  the  entire 
length  or  entire  duration  of  the  measurements  procedure. 

Discussion:  The  variance  <862>  in  electrical  terms  is  the  averaged  ac  power  over  a  specified  sam¬ 
ple  size  of  finite  dimension: 

c  ^ 

J0  <|S0|2>  dx  -  L  <862  > . 

It  is  related  to  the  power  spectrum  8S^{a)  of  the  random  process  d9 (x)  in  that  8S5#(a)  shows  how  the 
variance  is  distributed  over  wavenumber.  In  numerical  calculations  the  value  to  be  assigned  to  the 
space-averaged  variance  will  be  determined  in  this  report  by  integration  of  the  power  spectrum  over  a 
band  of  wavenumbers.  Because  L  is  finite,  the  power  contained  in  spatial  frequencies  whose 
wavelengths  are  much  larger  than  L  will  not  be  present  in  <S9[> .  Hence  the  allowable  band  of 
wavenumbers  must  be  finite.  Furthermore  there  is  an  additional  limit  on  the  thermal  energy  in  the  fiber 
over  the  length  L.  It  is  found  in  <  X2> : 

<X2>  =  <  fQJQ  &9(x)b9(x ')  dxdx'  > 

=  2  <b»2>  f0Lfo'  p(0  dtdx'. 

This  shows  that  lack  of  correlation  between  the  temperatuie  at  the  beginning  and  at  the  end  of  the  fiber 
reduces  the  thermal  energy  in  the  fiber. 

In  summary,  the  method  used  in  this  report  is  to  calculate  <b02>  over  a  finite  band  of  spatial 
wavenumbers  (or  temporal  frequencies)  and  then  to  evaluate  the  effect  of  correlation. 
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We  relate  these  results  to  those  of  Garrett  and  Munk  GM72  [31;  the  appropriate  formula  is  that 
of  their  footnote  on  page  237, 

12 


<802(Aa)>  —  f  Fr{a)  dot  —  f 

J  Aa  J  Aa  Qy 


/{(a)  tfa, 


in  which  the  units  of  Tv  (a)  are  (m2/cycle)/m.  Thus  the  mean-squared  phase  is 


<802(A a)>  =  C„2  2LLC  J 


</0 


dy 


FAa)  dot  (units:  rad2). 


Letting,  as  before,  2LLC  =  0(l)m2,  we  choose  a  fiber  of  length  L,  index  of  refraction  n ,  and  laser 
wavelength  X0  *  2 7t/A'0,  for  which  C#  is  known  (see  later  sections),  and  write 


<8<£2(A«)>  —  k0m 


1_  9w  J3L 

n  dO  L  dO 


V  <802fAa). 


In  these  formulas  the  integration  limits  enclose  a  band  of  wavenumbers  Aa  wide.  The  limits  can  be 
extended  to  plus  or  minus  infinity  if  a  window  W(x)  is  used  to  represent  the  effects  of  measurement 
restrictions  or  signal-processing  restrictions.  A  formal  discussion  of  windows  is  found  in  Appendix  B. 
For  the  purposes  of  this  report  we  choose  a  simpler  approach.  Let  W(x)  (units:  m-1)  be  a  spatial  win¬ 
dow  modeled  as  an  impulse  response  of  a  filter.  Then  the  Fourier  transform  of  this  filter  is 

W (a)  =  -t— -  f  lV(x)e~iaxdx  (units:  none). 

ITT 


Thus  the  band-limited  variance  of  8</>  is 


DO 

dO 

l-oo 

dy 

F7 (a ) |  fF (a ) | 2  da  (units:  rad2). 


In  a  similar  way,  if  w(r)  is  a  temporal  window,  then 

W  (tu)  =  — —  f  w(t)e~,<otdt 

2tt  j 

and 


<b<t>j>  =  2  C„2r  7T(  / 


d9 


dy 


F((a>)  |  W(u>)\2  dco, 


where  F-(cu)  is  the  temporal  power  spectrum  of  internal  waves  and  I  W{u))\2  is  the  filter  transfer  func¬ 
tion. 


The  procedure  of  analysis  summarized  above  depends  very  strongly  on  the  existence  of  measured 
power  spectra  of  internal  waves  in  the  ocean.  The  precision  of  these  measurements  and  their  ultimate 
significance  under  the  conditions  of  local  space  and  time  make  them  only  approximate.  In  most  cases 
good  approximations  are  available.  These  will  be  used  in  later  sections  of  this  report  to  made  numerical 
estimates  of  environmental  noise,  which  will  ultimately  be  used  to  determine  a  signal-to-noise  figure  of 
merit  S/N  of  the  fiber-optic  sensor. 

The  mathematical  form  of  S/N  points  the  way  toward  the  appropriate  form  to  be  given  to  the 
environmental  noise  effect.  We  now  turn  to  an  investigation  of  the  required  general  expression  for 
S/N  of  a  fiber-optic  hydrophone. 

SIGNAL/NOISE  RATIO  OF  FIBER-OPTIC  HYDROPHONES 
Photodetector  Current 

In  a  laser  heterodyne  detection  system  the  superposition  on  the  surface  of  the  photodetector  of 
the  reference  beam,  with  intensity  /r,  and  the  signal  beam,  with  intensity  /s,  can  be  modeled  on  the 
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theory  of  partially  coherent  light  |4).  In  this  theory  the  total  intensity  of  two  superimposed  beams 
which  differ  in  phase  by  A <t>  is 

In  =  /,  +  Is  +  2|y rs\y[TJs  cosier,,  (r)  -  A<£]  (units:  V-A/m2), 
in  which  yrs  is  the  complex  degree  of  mutual  coherence,  lyrs|  is  its  modulus,  and 

ars  (r)  =  argy„  + 

where  f  is  the  mean  frequency  of  the  laser  and  r  is  spatial  path  difference.  In  terms  of  the  laser-beam 
wavenumber  A,  one  can  write  the  phase  difference  as  a  path  difference  Ax, 

A<£  =  k  Ax  (units:  rad). 

If  the  laser  field  is  represented  by  the  analytic  signal  V{t ),  then  by  definition 

=  rr5(r)  =  rr5(r) 

y “  T  VVo)  Vivo)  Va  VA 

and 


r„(r)  =  <  Vr(t  +  T )  v;u)>  (units:  V-A/m2). 


in  which  <  >  signifies  statistical  average  in  time,  V rr  is  the  autocorrelation  of  the  random  variable 
V(t),  Vrs  is  the  crosscorrelation  of  the  random  variables  Vr  and  and  r  is  the  delay  time.  Since  \yrs  I 
cos  [a,, (r)  +  A </>]  is  a  real  number,  Reyrs,  it  is  a  measurable  quantity,  namely, 


Rey« 


I  rs  zJiZJl 

2  Va  Va 


The  total  intensity  !n  thus  is  a  sum  of  coherent  and  incoherent  parts.  It  is  justifiable  to  assume  that 
over  sufficiently  long  time  the  random  part  is  uniformly  incident  over  the  area  of  illumination  of  the 
photodetector.  By  integration  over  this  area  the  total  light  power  incident  on  the  collecting  surface 
becomes  J 

Pn  =  Pr  4-  Ps  +  2 \yrs\JPrPs  COs[an(t)  —  A<£]  (units:  V-A). 


This  light  power  generates  a  photodetector  current  of  time-average  value 

</>  =  D<Prs>  (units:  A), 

in  which  D  is  the  photodetector  conversion  factor  (units:  V”1  or  A/VA).  To  emphasize  partial 
coherency,  it  is  useful  to  write  this  current  in  the  form 

</>  -  |yj  D[Pr  +  Ps  +  2  cos(aw  -  A<^>)]  +  (1  -  \yn\)  D(Pr  +  Ps). 

The  first  term  represents  the  coherent  superposition  of  two  beams  of  powers  |yr5|  Pr  and  |yrJ  P5,  while 
the  second  term  represents  incoherent  superposition  of  two  beams  of  powers  (1  —  |y„|)  Pr  and 
(1  -  | y r5 1 )  Ps .  This  distinction  between  coherent  and  incoherent  parts  of  the  photodetector  current  will 
be  used  to  advantage  in  further  developments. 


The  phase  A <f>  is  the  sum  of  two  terms:  a  nonacoustic  phase  A0O  and  an  acoustically  induced 
phase  In  the  absence  of  the  latter  the  average  photodetector  current  can  be  written  just  as  shown 

above  with  A <t>  replaced  by  A $0-  It  is  then  essentially  a  sum  of  dc  terms,  Iph  (units:  A).  In  the  pres¬ 
ence  of  an  acoustic  signal  which  varies  sinusoidally  in  time,  the  cosine  term  becomes 

cos(a„  -  A0O—  A<£4)  =  cos(ar5  -  A$0)cosA<£4  +  sin(ars  -  A </>0)  sinA^. 

If  we  assume  that  sin  A<£4  =  A 04,  then  there  exists  a  time-varying  photodetector  component  current, 

<iAc>  =  2D\yrs\  JPrPs  sin (a„(r )  -  A<£0]  A <t>A  U). 

The  nonacoustic  phase  A<£0  may  contain  a  random  component  of  time  A0O(;)  caused  by  environmental 
disturbances.  In  this  case,  when  there  is  no  useful  acoustic  signal. 
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cos[a„  -  A0O“  A</>0(r)]  =  cos(a„  -  Ad>o)cosA0o(/)  +  sin(ars  -  A0o)sin  A  <A0  (/)• 

When  the  fixed  phase  ars  -  A <t>o  is  adjusted  to  be  n/2 ,  as  is  customary,  there  remains  a  random  com¬ 
ponent  of  photodetector  current 

4*0  =  2Z)  lyj  sin  [A0o(/)I. 

The  probability  distribution  of  is  determined  by  the  sine  of  the  probability  distribution  of  A<£o(/). 
In  later  sections  of  this  report  we  shall  consider  A<£0(>)  to  be  caused  by  random  temperature  fluctua¬ 
tions,  and  we  shall  derive  numerical  estimates  of  its  magnitude  in  typical  ocean  environments. 

Coherent  and  Incoherent  dc  Photodetector  Currents 

The  time-averaged  dc  photodetector  currents  contribute  noise  in  the  detection  circuits.  As  noted 
above  these  currents  are  the  sum  of  coherent  and  incoherent  terms.  The  incoherent  terms  consist  of 
three  contributions: 

•  scattering  of  laser  light  in  the  fiber  b,  Bragg  reflection  from  thermally  induced  acoustic 
plane-wave  trains  (Brillouin  scattering); 

•  scattering  from  mass-density  fluctuations  caused  by  random  motion  of  the  molecules  of  glass 
in  the  fiber  whenever  the  wavelength  of  the  incident  laser  beam  is  much  smaller  than  the 
size  of  the  molecules  (Rayleigh  scattering);  and 

•  scattering  from  spherical  (or  irregular)  inhomogenities  where  the  characteristic  dimension  is 
larger  than,  or  comparable  with,  a  wavelength  of  the  incident  light  (Mie  scattering). 

Thus,  the  incoherent  part  of  the  dc  current  may  be  explicitly  written  in  the  form 
D(  1  -  lyjHP,  +  Ps) 

=  D  (PrH  +  PrK  +  Pr\f  +  I\„  +  PsR  +  PsM)  (units:  A). 

where  the  subscripts  B,  R,  and  M  refer  to  Brillouin,  Rayleigh,  and  Mie  scattering,  respectively.  The 
coherent  part  of  the  dc  current  has  already  been  formalized  above: 

(  ^ ph  ^  coh  \D[Pr  +  Ps  +  2VKK  cos(a„  4-  Ac£0)l  (units:  A). 

Here  all  terms  are  non-time-varying  ones. 


Noise  in  the  Photodetector  and  Associated  Preamplifier  Circuitry 

Shot  Noise 

Shot  noise  power  is  determined  from  the  dc  current  lph  of  the  photodetector  explicitly  formulated 
above  and  two  other  currents:  a  dc  current  due  to  background  light  (/#*)  and  the  dc  dark  current  which 
enters  the  detector  circuit  (//,).  If  the  photodetector  has  an  internal  current  gain  G  the  total  shot  noise 
will  be  given  by 

2<'((/,,Ai,h  +  l  lph) incoh  +  4a,  +  //>]  <G2> 

in  which  e  is  the  electronic  charge  and  <<j2>  =  <G>2  F{G),  where  F(G)  accounts  for  an  additional 
multiplication  (or  gain)  due  to  the  nonlinear  properties  of  gain  development  inside  the  photodetector. 
The  value  of  HG)  is  0.5  to  3.0,  depending  on  the  type  of  detector  used  15], 


units: 


Hz 
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Thermal  Noise 

The  load  resistor  RL  of  the  detector  generates  Johnson  noise  Wj  when  a  thermally  induced 
current  ir  flows  through  it: 

Wj  =  <//>  Rl. 

This  noise  power  in  a  band  8/  units  wide  at  absolute  temperature  0  is 

(units:  V-A), 

where  k  =  1.38  x  10  23  J/K  (Boltzmann's  constant).  The  noise  density  is  thus 

A2' 


<Jt>  _  4*0 
8/ 


Ri 


units: 


Hz 


This  noise  component  is  additive  to  the  shot  noise  formulated  above. 
Laser  Jitter  Noise 


Laser  beams  have  associated  vith  their  steady  intensity  /  a  corresponding  fluctuating  intensity  8/. 
This  is  the  laser  jitter.  It  gem  s  in  the  photodetector  a  noise  current  ij(t).  If  a  data  record  T 
seconds  long  of  this  current  is  Courier  analyzed  it  will  have  a  power  spectral  density  \/j(co)\2  and  a 
power  spectral  density  factor,  defined  here  to  be  the  quantity 

in  which  Ij(oj)  (units:  As)  is  the  Fourier  transform  of  if(t)  over  the  duration  of  the  finite  sampling 
time.  The  choice  of  Fourier  transform  is  explained  in  Appendix  C.  The  mean-squared  value  of  //(/)  is 
then  obtained  by  integration  over  the  band  of  frequencies  spanned  by  Ij(u>)  in  the  power  spectral  den¬ 
sity  factor: 

<ij( Aa>)>  =  f  \fj(to)\2  dw  (units:  A2). 

•r  Auj  / 

A  convenient  way  of  deriving  the  power  spectral  density  factor  is  through  use  of  the  autocorrelation 
Y j(t)  of  ij(t): 

~\lj(w)V=  -7-  C  eim\(r )  dr. 

T  ATT 

where 

V(r)  =  y  J*r  ij(t)  ij  (t  +  r)dt  (units:  A2). 

Typical  forms  of  Y 'j(t)  are  discussed  in  a  later  section  of  this  report. 

Environmental  Noise 


Assume  that  the  environment  generates  a  random  thermal  disturbance  and  a  random  mechanical 
disturbance  of  the  fiber.  Let  one  of  these  disturbances  be  8e(r,f),  meaning  it  is  a  spatial  and  temporal 
fluctuation  of  an  environmental  parameter  (temperature  or  mechanical  stress).  The  power  spectral  den¬ 
sity  factor  of  this  random  function  over  finite  volume  V  and  finite  time  Fis  then  defined  to  be 


(2tt)a 


VT 


in  which  the  quantity 


00 

E( K,  w )  =  ff  Mr.  t 


dx 

(2tt)3 


7 


w 
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is  the  Fourier  transform  of  be  (r,  /)  (Appendix  C).  The  mean-square  value  of  be  over  these  intervals  is 

<8«  2>  =  //  (2tt)4  </3K  dto. 

A  convenient  way  of  deriving  the  power  spectral  density  of  be  is  through  use  of  the  autocorrelation 

Y*€  (d,  r): 

YS{(d,T)e'“'-'Kd^r</d=  (2tt)4  ^£(K^-)-1-, 

where 

i  00 

Yfi((d,  r)  -  -j~J J  8e  (r,  /)8e  (r  +  d,  /  +  t)  dl  di d. 


The  power  spectral  density  per  unit  of  bandwidth  can  be  derived  from  these  formulas  by  integra¬ 
tion  over  wavenumber: 

~\E(to)\2  =  f  (2tt)3  *£-^:/)*2  d3 K. 

By  a  linear  conversion  factor  the  units  of  this  parameter  can  be  made  A2-s.  In  most  applications  the 
spatial  part  of  the  spectrum  is  considered  essentially  one-dimensional.  Then  the  power  spectral  density 
factor  of  be  (x  t)  becomes 

(2n)2  1  E(a,(o)  l2 

LT 

where 

E(a,u) )  -  f  f  8«U  t)eiu,,~iax-^-~. 

J  J  ATT  ITT 

The  mean-square  value  of  be(x,  /)  is  obtained  from  this  by  integration: 

<8e2>  =  ff  (2tt)2  dadu). 

Also  one  can  form  the  power  spectral  density  factor  per  unit  of  frequency  bandwidth, 

~r  |£(w) I2  =  f  (2t r)i^«liirfa. 

Again  by  a  linear  conversion  factor  the  units  of  this  parameter  can  be  made  A2  s. 

In  all  cases  the  quantity  V <8€2>  is  the  desired  representation  of  the  environmental  disturbance, 
averaged  over  space  and  time.  The  root-mean-square  amplitude  of  fluctuation  in  the  phase  of  the  laser 
signal  beam  is  then  desired  to  be 

A<*>rms  =  A(-J <8e2> ,  A€  =  const. 

However,  since  the  goal  of  this  analysis  is  to  estimate  the  noise  current  in  the  photodetector  caused  by 
environmental  disturbances,  the  quantity  actually  calculated  is 

<i2>  =  Bf  < be2> ,  Bf  =  const, 

in  units  of  A2.  For  purposes  of  relating  environmental  noise  to  other  noises  described  above,  it  is  con¬ 
venient  to  represent  it  in  the  form  of  noise  per  unit  bandwidth: 

</f2>  x 
— - —  Ce  <  be2( (o)>,  Ct  =  const. 
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\ 

f  V 


Examples  of  the  application  of  these  formulas  to  environmental  disturbances  which  affect  the 
transmission  of  light  in  fibers  are  given  below.  We  consider  here  the  effect  of  these  disturbances  on 
the  phase  of  the  laser  light. 

Let  L  be  the  length  of  a  fiber  which  is  acting  as  a  transmission  path  for  a  laser  beam.  As  the  light 
travels  from  the  origin  to  L  its  phase  is  altered  by  an  amount 

0  =  fQ  ko(x)n(x)dx  (units:  rad), 


where 


k0ix)  =  background  wavenumber  =  27 r/A0(x)  (units:  m  *) 


nix)  =  index  of  refraction  (units:  none). 

A  simple  case  occurs  when  k0n  is  a  constant  along  the  transmission  path.  Then 

0  =  konL. 

Now  assume  that  the  fiber  is  acted  upon  by  a  random  environmental  disturbance,  say  mechanical  stress, 
which  is  represented,  for  convenience  in  this  report,  as  a  scalar  pressure  8SP  (units:  N-m-2)  normal  to 
the  surface  of  the  fiber,  or  as  a  scalar  torsion  6 Sr.  The  change  in  phase  to  first  order  in  the  change  in 
stress  is 

80  * =  Si  * + *°ms + k°nL  Os 8s  (unks:  m_,)’ 

in  which  the  asterisk  indicates  "per  unit  length”  For  scalar  pressure  changes  the  last  term  in  the 
parentheses  can  be  written  in  terms  of  the  bulk  compressibility, 

s  ~  Ar  (units:  m2-N~!). 
oop  V 

If  AL  represents  the  change  in  one  dimension,  then  for  a  three-dimensional  isotropic  body 


jr1-- 


1  AL 
hSp  L 


Replacing  dL/dSp  by  A L/dSp,  and  assuming  dk(Jd8Sp  ==  0,  one  finds 


=  MI  =  iii  _  E_ 


,  rad-  m2 

k0n  units:— - 

N  •  m 


Following  the  procedure  explained  earlier  in  the  section  "Procedure  of  Analysis"  the  average  value  of 
the  fluctuation  in  60  is 

f><t>LTms  =  V <80/; >  =  cj <dSp>  2 LLC 

( 1  Bn  K-'\  ,  n^Trr- 


=  \  n  S - 3~ 'j  k$n  v  <8SP2>  (units:  rad), 

assuming  that  2 LLC  =  0(1)  m2.  A  similar  formula  can  be  constructed  for  changes  in  scalar  torsional 
stresses.  We  assume  that  a  constant  twisting  moment  M  induces  a  constant  angle  of  twist  0  in  a  fiber 
of  circular  cross  section  with  polar  moment  of  inertia  /o, 

M  -  Gl0j-  (units:  N-m), 

in  which  G  is  the  modulus  of  rigidity  (units:  N-m'2),  and  the  units  of  70  are  m4.  A  change  in  M 
induces  a  change  in  0: 


f 
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MU  =  GI0^-. 

We  assume  a  linear  relation  exists  between  a  change  in  fiber  length  and  a  change  in  angle  of  twist.  By 
definition, 

a  =  '~r  (units:  m). 

A0 

Hence, 

A  L  _  a 
LkM  ~  GIq  ’ 


Thus,  for  random  twisting  moments  8M, 


Wums  =  V<Hl>  =  --|£  +  M  V <8A/2>  2LLc. 

n  a  M  UIq 


This  is  a  one-dimensional  statement.  In  three  dimensions  we  repeat  the  above  analysis  to  obtain  the 
effects  of  random  torsional  stress  SST  : 

v<8<*2>  =  M  V <&s) >  2llc, 

n  OJj-  3 


A'ow  V  <8A'72  >  2Z-Lr, 


in  which 


In  both  cases  the  symbol  8  is  used  to  indicate  fluctuations.  The  constant  Cs  is  directly  related  to  special 
constants  in  the  theory  of  photoelasticity,  values  of  which  appear  as  C  in  Part  II  of  this  report,  NRL 
Report  8425. 


Laser  Jitter  and  Environmentally  Induced  Noise  for  Gaussian  Correlated  Fields 

In  the  section  "Environmental  Noise"  it  was  found  that  the  noise  current  in  the  photodetector 
current  can  be  obtained  from  the  mean-square  environmental  disturbance  <5e2>  averaged  over  space 
and  time.  This  quantity  in  turn  can  be  calculated  from  a  knowledge  of  the  power  spectral  density 
|£(K,  oj)|2  of  €  (K,  /).  Again,  as  noted  earlier,  a  convenient  way  of  calculating  spectra  is  through  the 
use  of  autocorrelation  functions.  We  treat  now  several  cases. 

Power  Spectral  Density  of  Mechanical  Stress  and  Associated  Noise 

Let  the  environmental  disturbance  be  a  random  mechanical  stress  which  is  represented  in  reduced 
form  as  a  scalar  fiSTr,  t).  The  autocorrelation  function  over  all  space  and  time  is 

\  S(L  t)  =  lim  f  f  hS (r,  /)8S(r  +  1,  t  4*  r)dvdt  . 

5  k— oo  VT J  v  J  / 

1-00 

In  all  applications  to  real  phenomena  both  V  and  T are  finite.  The  mean-square  fluctuation,  or  variance , 
is  the  value  over  the  whole  Land  the  whole  T : 

<8S2>  -  Y,(0.  0)  =  f  f  8S2(r,  t)  dvdt  (units:  N2  m"4). 

VT  J  T 

Because  of  instrumental  and  experimental  limitations,  not  all  space  and  not  all  time  can  be  measured. 
There  is  an  effective  band  of  wavenumbers  and  an  effective  band  of  frequencies  that  are  resolvable  under 
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the  conditions  of  test.  Thus  the  full  variance  <852>  is  not  physically  reportable.  To  deal  with  a  finite 
spatial  and  temporal  spectrum,  we  utilize  the  full  autocorrelation  function.  Examples  are  now 
developed. 

Assume  that  the  scalar  stress  field  is  isotropic  in  space,  so  that  its  autocorrelation  is  dependent  on 
distance  /  between  r  and  r',  and  not  on  orientation  of  the  vector  separation  I.  Further  assume  a  simple 
Gaussian  form  for  both  space  and  time  dependence: 


1 

/  V 

1  (  T  2 

2 

Km 

2  T,m  j 

in  which  henceforth  we  shall  take  the  units  of  <852>  to  be  the  noise-current  equivalent,  A2.  Here  lim 
is  the  correlation  length  of  the  random  stress  field  and  rcm  is  its  correlation  time.  As  noted  earlier,  the 
power  spectral  density  factor  is  given  by 

lid  II  r  P 


(2tt)4  -§-5^aj)^  =  f  f  circle  2  c  2*T<WI  e"*  Vwr, 

VT  (h rr 

in  which  /a  is  a  vector  wavenumber.  The  integration  over  r  yields  \Fhr  rcm  exp[-  The 

integration  over  1  is  found  by  noting  that  =  dLxdLydLz  and  that,  because  of  isotropy  of  the  stress 
field  ,  dH  —  (\/277/cm)  V/v,  where  dlx  is  nondimensional.  Thus, 

(2tt)4 =  <852>  exp  --y I/a  pt|/2m  -  -yw2r2„,  (units:  A2m3s). 

By  integration  over  a  volume  of  effective  wavenumber  A fx  we  obtain  the  noise  power  spectral  density 
factor  per  unit  bandwidth  of  frequency  . 


27T  I  \  1 2  /jmrci 

— Issw)l  — sp— 


J*  1  1  ■>  ■> 

,  exp -  —  Im  n\lL  dln  exp  wJT,?m  (units:  A: •  s). 
V  2  2 


in  which 


|8S(W)P  =  f  (units:  A2  .s2) 


Finally,  by  integration  over  an  effective  band  of  frequencies  Aw,  one  gets  the  mean-squared  band- 
limited  variance: 


<852(A/a.Aw)>  - 


<8 S2>l?m  r„ 


J>  fs  (/V cxp  -\w  iim~  7oj’T«2^ 


- - - (27r)2  Ji/u  |  2  J 

=  <852>  F4(A/a.  Aw)  (units:  A2). 

Here  F4  (A /a.  Aw)  <  1  when  A /a  and  Aw  are  finite,  and  /^(A/a.  Aw)  =  1  when  A/a  and  Aw  are  infinite. 
The  subscript  4  denotes  a  four-coordinate  system  (3-space  and  time). 


When  the  spatial  dependency  of  85  is  one-dimensional,  the  power  spectral  density  factor  is 


<852>  .  Km 

~Th^r cmTim  exp  “  2 


(units:  A2*nvs). 


Thus  the  power  spectral  density  factor  per  unit  bandwidth  of  frequency  is 

11 


f 
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Band-limited  Variance  of  Laser  Jitter  Noise 

Assume  that  laser  jitter  is  predominantly  a  temporal  phenomenon.  Then,  using  the  equations 
derived  above  for  Gaussian  autocorrelation  of  the  time  record  of  jitter  current  in  the  photomultiplier, 
we  estimate  the  band-limited  variance  of  this  record  to  be 

<8J2>  t  — 

<8/y  ( Aw)>  =  - r-—’-  J  f  e  2  da)  (units:  A2), 

v  2t t  j 

where  tcJ  is  the  correlation  time  of  the  record  and  Aw  is  the  allowed  bandwidth  of  the  measurement 
system. 

Band-limited  Variance  of  Fluctuating  Temperature  Noise 

Assume  that  the  random  temperature  field  is  one-dimensional  in  space  and  possibly  fluctuates  in 
time  and  that  both  spatial  and  temporal  autocorrelation  are  Gaussian.  Then  the  band-limited  variance 
of  the  temperature  fluctuation  is 

<802(Ac*,  Aw)>  —  — — —  lc9rc#  f  dco  e  2  f  dae  2  , 

ITT  **  Aw  A« 

where  lcii  and  tcB  are  the  correlation  lengths  of  the  spatial  record  and  the  temporal  record  respectively. 
By  use  of  an  appropriate  conversion  factor  this  becomes  also  the  band-limited  variance  of  the  noise 
current  generated  by  environmental  noise. 

The  use  of  Gaussian  autocorrelations  for  the  random  disturbing  field  may  not  be  a  best  choice.  In 
application  to  ocean  environments  a  sin  x/x  autocorrelation  has  the  possible  advantage  of  greater  fidel¬ 
ity  to  the  random  phenomena.  A  derivation  of  power  spectral  density  factors  in  the  sin  x/x  case  is 
presented  in  Appendix  D. 

Noise  Sources  in  the  Preamplifier  [5] 

A  preamplifier  can  be  schematically  represented  as  a  noiseless  generator  with  gain  A  (w)  and  input 
admittance  YA  =  GA  +  j(o C A .  It  has  two  noise  sources,  a  voltage  noise  source  vA  in  series  with  the 
input  terminals  and  a  current  noise  source  in  parallel  with  the  input  terminals.  The  latter  consists  of  a 
portion  iA  which  is  uncorrelated  with  the  noise  voltage  source  and  a  second  portion  which  is  correlated 
with  this  source  through  an  admittance  Y,  =  ReK(  +  j  Im  Yc .  We  neglect  this  second  portion.  The 
spectral  density  of  noise  from  the  voltage  source  is  <vA>/8f\n  units  of  V2/Hz,  and  that  of  the  current 
source  is  in  units  of  A2/Hz. 


Noise  Voltage  at  the  Output  of  the  Preamplifier 

The  noise  currents  in  the  photodetector  and  preamplifier  circuits  are  additive.  Because  of 
different  methods  of  derivation,  some  of  these  noise  densities  are  on  a  unit  bandwidth  6/ basis  (in 
cycles/s),  while  the  remainder  are  on  a  unit  angular  "adian  8w  basis  (in  rad/s).  It  will  be  convenient  in 
this  report  to  keep  them  in  separate  groups.  Thus,  the  noises  NhJ  and  Nhio  are: 

On  a  8/  basis:  Nh  f  =  shot  noise  +  amplifier  noise  +  Johnson  noise. 

On  a  8w  basis:  N ^  —  laser  jitter  noise  4-  environmental  temperature  noise 

+  environmental  mechanical  stress  noise. 
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The  total  noise  in  units  of  V2  at  the  output  of  the  preamplifier  is  the  sum  of  two 
Ntf  =  fs/  2 «-[</*>«*  +  (4J,ncoh  +  Ibk  +  h\  <G>2F(G)  + 


terms: 


and 


<v}> 

8/ 


Ik,  +  yc \2  + 


AkT 


Rl 


\Ai o>)l2 

IK/+  Kj2 


df 


8/ 

(units:  V2) 


N 


but 


- ^  [i*_  |SS/,„)|. w„(„. +  if  ;  faulty. «>i1.  d, 


U(co)|2  , 

X  - TT 

Y,  4-  V,  2 


T 

(units:  V2), 


where 


ISS/M2  is  the  power  spectrum  of  laser  jitter 
|§SM(a,  w)|2  is  the  power  spectrum  of  random  temperature  effects 
co)\2  is  the  power  spectrum  of  random  mechanical  effects 
and 

K0  is  the  volume  of  space 

If  we  assume  Gaussian  autocorrelations  of  the  random  disturbances,  then 


(units:  A2-s2), 
(units:  A2  •  m2  •  s2), 
(units:  A2  -  m6*s2), 

(units:  m3). 


r  <8J2>  tcj 
Jaw  V2t t 


<S02>Il„tc 

2tt 


L 


da 


<852> 

(In)2 


li  St  T  CM  J, 


v  d'fiex p[-—  •  nUcM  ~ 


\A  (ou)|2 

Ik,  +  Kfl2 


d<o. 


Thus,  in  order  to  calculate  from  assumed  Gaussian  autocorrelations,  we  are  required  to  know  the 
total  variances  of  the  random  laser  jitter,  temperature,  and  mechanical  stress  field,  as  well  as  their  tern- 
poral  and  spatial  correlation  lengths. 


Signal  Current  and  Voltage  Sensitivity  of  the  Hydrophone 

In  the  section  "Photodetector  Current"  it  was  seen  that  there  is  a  lime-varying  component  of  pho¬ 
todetector  current,  which  we  call  here  the  signal, 

^signal  2D\yrs  \  yfp—ps  sin  [ars  (r )  -  A<£0]  A 4>AU), 

where  A is  the  phase  generated  by  the  acoustic  signal,  considered  to  be  of  small  enough  magni¬ 
tude  such  that  sin  A <f> A  =  A <t>A.  An  oft-used  design  feature  is  to  experimentally  adjust  ar5(r)  -  A<£0  to 
be  an  odd  multiple  of  W2.  This  procedure  maximizes  the  signal.  Assume  that  this  is  done. 

Now  assume  that  A</>  4  (i)  =  Ad)^  sino>5 1.  Then  the  rms  ac  component  of  photodetector  current  is 

( ^signal) rms  ~  2D\yrs\JPTPM^Ahms- 

The  change  in  phase  is  linearly  related  to  the  acoustic  pressure  amplitude  pA.  Hence  (from  earlier  sec¬ 
tions  in  this  report)  one  has 
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(/ 


signal' rms 


2D\yrs\JPTs 


-L 

n  bp 


*0 nL  I pA  lrms. 


The  result  is  characterized  by  the  conversion  of  phase  modulation  into  amplitude  modulation.  Actually 
the  signal  current  is  a  statistical  average,  since  the  detection  scheme  is  assumed  only  partially  coherent. 
We  therefore  write  it  as  an  average,  <  ( /‘signal rms- >■ 


In  photodetector  procedures  the  signal  current  appears  across  the  photodetector  load  resistor  Rt  . 
where  it  develops  the  signal  voltage.  Hence  the  receiving  voltage  sensitivity  M  is 


M  = 


<(WJg>^  =  2DRl  )yj  JiT]? 


IP  ir 


n  dp 


knnL 


(units:  V  •  m2 


N  '). 


This  is  the  sensitivity  without  a  preamplifier.  The  effect  of  the  preamplifier  is  discussed  next. 

Signal  Intensity  at  the  Output  of  the  Preamplifier 

In  formulating  the  S/N  ratio  of  the  hydrophone  we  have  elected  to  write  both  S  and  N  in  units  of 
V2  which  appears  across  the  load  resistor.  As  noted  earlier  the  photodetector  itself  amplifies  the  current 
by  amount  <G>.  The  preamplifier  also  amplifies  the  signal  by  amount  A  (a>).  Hence  the  signal  at  the 
output  of  the  preamplifier  has  an  intensity 

5=  <(/5,g„a,)rms>3  <^>2  J^Tjl  (“nits:  V2). 


where  Y/  and  YA  are  component  electrical  admittances  of  the  photodetector  amplifier  circuit  (Appendix 
E).  It  is  to  be  noted  that  the  signal  is  single  frequency.  If  the  signal  has  the  character  of  banded  noise 
it  can  be  reported  on  a  per  bandwidth  basis  <i2>/hf.  Then  the  signal  intensity  appears  in  the  form 

[a  («j)P 


-L 


<i}> 


<G>' 


*/  ""  I  Y,  +  Ya  \2 


df 


(units:  V2). 


Explicit  Form  of  S/N  Ratio  of  the  Fiber-Optic  Hydrophone  at  the  Ouput  of  the  Preamplifier 

For  the  convenience  of  the  reader  we  repeat  here  the  explicit  form  of  the  Sf  N  ratio  of  a  fiber¬ 
optic  hydrophone  for  the  general  case  of  a  banded  noise  signal: 


N 


s. 


<4  > 

4/  8/ 


<G>‘ 


\j4  (to)  |2 

\Y,+  Ya  I2 


df 


ff  2c[(/p„)C0h  +  (4/,)incoh  +  4a  +  41  <G>2  F(G)  + 


<Ja> 

8/ 


+  U7+  Yc\2+jf 

oj  K  L 


I A  (<u)  |2  ... 

I y,yaV  3 


+  J-y-  \8Sj(o,)\2  +  f  ^-lSSs„(a,(o)l2da 
2n_  C  (2n)3 

T  J  V« 


I «)Vd^ 


\_A  (co)l2 
I  Y,  +  Y4 12 


All  symbols  have  been  defined  in  previous  sections  of  this  report. 
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INSTRUMENT  WINDOWS  AND  THE  NYQUIST  LIMIT 

The  random  environmental  disturbances  which  affect  the  transmission  of  light  through  a  fiber 
have  autocorrelation  functions  which  die  away  naturally  with  increasing  spatial  and  temporal  lags. 
Hence,  the  power  spectral  densities  have  continuous  spectral  spreads  from  — 00  to  +°°.  The  variances 
of  the  random  disturbances  can  be  calculated  by  integration  of  the  spectrum  over  time-infinite  limits. 
However,  because  of  instrument  limits  in  time  and  space,  and  because  of  the  conditions  of  field  experi¬ 
ments,  only  a  finite  band  of  frequencies  and  a  finite  band  of  wavenumbers  are  actually  measured.  Thus 
only  a  portion  of  each  variance  of  a  random  disturbing  field  appears  in  the  averaging  process.  An 
important  problem  of  analysis  centers  on  determination  of  the  effective  frequency  and  wavenumber 
bands.  We  discuss  this  next. 

Two  choices  of  Fourier  space  are  in  common  use.  In  the  first  the  Fourier  pairs  are  ( k ,  x)  and 
(w.  /),  and  in  the  second  the  pairs  are  (A-1,  x )  and  (/,  /).  To  illustrate,  choose  a  function  0(x,  t ) 
which  is  one  dimensional  in  space  and  which  varies  with  time.  Give  it  the  dimension  of  temperature, 
°C.  Then 

S„(k.  •>-  ff  (Hx.  t)em,-,kx  ~~~rx  (units:  °Cms) 

J  J  (2tt  )  * 

and 

SH  (A-1,  /)  =  ff  0(x,  ()e,2nf,-,2irK  '*  dtdx  (units:  °C*ms). 

Let  us  use  (A-1,  /)  space  to  illustrate  a  first  instrument  limit  on  the  spectrum.  Consider  one  coordinate 
at  a  time: 

S„(/)-J  U(t)e'2”"dt. 

Suppose  0(t)  is  sampled  at  intervals  A,;  then  the  modulated  signal  is 

0,  (/)  =  ti(t)  £  6(/  -  nA,) 

00 

and 

S,(J  )  =  ~~  £  S„ 

-00 

Thus  the  sampled  time  signal  has  a  transform  with  period  1/A,;  that  is,  the  Fourier  transform  repeats 
itself  indefinitely.  If  the  frequency  spread  of  the  continuous  S#(/)  overlaps  the  period,  then  frequency 

components  S{f)  for  /  >  appear  in  adjacent  spectra.  Thus  the  limiting /for  complete  recovery 

of  0(t)  is 

fs  =  “tt-*  (units:  cycle/ s). 

ZA, 

This  is  the  Nyquist  limiting  frequency.  In  a  similar  way,  for  the  spatial  function  0(x)  sampled  at  inter¬ 
vals  1/2  A,,  the  spectrum  is  cut  off  at 

A/J1  =  -r-j—  (units:  cycle/m). 

2A4 

The  Nyquist  limits  in  (k  —  o>)  space  are 

1  7T 

ut\  ”  2 nfs  =27 r  x  — —  —  —  (units:  ra<Vs) 

2A,  A, 
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and 

kN  -  2irX*1  =2 ir  x  (units:  rad/m). 

If  the  measurement  instrument  samples  time  and  space  continuously,  the  Nyquist  limits  are  infinite  in 
all  cases. 

The  Nyquist  limit  refers  to  sampling  the  spatial  and  temporal  record  of  the  environmental  field 
either  during  the  experiment  or  during  the  signal  processing.  We  discuss  next  other  instrument  and 
experiment  windows. 


a  —  <o  Power  Spectrum  and  Cutoffs 

It  was  noted  earlier  in  this  report  that  Rayleigh’s  theorem  allows  us  to  write  the  total  energy  in  a 
system  over  real  coordinates  as  the  integral  of  the  absolute  value  squared  of  the  spectrum  of  the  system 
over  all  spectral  components.  Applying  this  to  the  selected  case  of  a  fluctuating  temperature  field,  we 
note  that  the  total  energy  so  calculated  is  the  variance  of  the  fluctuating  field.  However,  because  of  the 
nature  of  the  measurement  procedure  only  a  portion  of  this  variance  is  actually  picked  up  by  the 
receiver.  This  portion  is  enclosed  within  the  bounds  of  the  various  cutoff  frequencies  and  cutoff 
wavenumbers.  Figure  1  shows  this  schematically. 


As  noted  earlier,  the  band-limited  variance  is  given  by 

<f>92( Aa,  Aco)>  =  f  f  | (a,  a>)|2  dadb), 

J  Aa  J  Aw 


where 

oo 

5 SM  (a.  oj)  =  — T-  JT  80U  /  dtdx. 

We  consider  the  cutoff  frequencies  of  the  fluctuating  temperature  field  80 (x,  t)  first.  Assume  that  the 
sensor  is  a  spatial  point.  The  field  80(x0.  t )  can  then  be  decomposed  into  its  ca-spectral  components. 
Let  the  receiver  have  a  finite  response  time  A  Tr.  Since  80(i)  is  a  random  process,  some  portion  of  it 
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will  generally  have  temporal  excursions  which  are  shorter  in  duration  than  A  Tr.  The  sensor  will  not 
respond  to  them.  Thus,  the  spectrum  of  80  will  be  cut  off  at  a  frequency  o>max  such  that 


Frequencies  higher  than  <omax  will  not  effectively  contribute  energy  to  the  receiver.  Of  the  remaining 
band  0  to  wmiiX,  a  second  portion  will  be  cut  off  by  the  integration  time  A  T{  allowed  for  the  receiver  to 
measure  the  field.  This  time  is  purposely  limited  to  permit  a  high  rate  of  data  collection  and  to  prevent 
the  random  field  89  from  drifting  into  nonstationarity.  The  cutoff  at  the  low  frequency  end  is  then 
estimated  to  be 


Between  the  two  cutoffs  one  then  finds  the  effective  frequency  band  of  the  to  spectrum. 

Aw  =  a>max  -  (Omin  =  aV  ~  aV’  A7^<A7'- 

In  towed  systems  the  upper  frequency  cutoff  of  the  internal  wave  spectrum  may  be  controlled  by  the 
speed  of  tow.  This  is  discussed  later. 

We  next  turn  to  the  determination  of  the  cutoff  wavenumbers.  Assume  that  the  receiver  is  a 
linear  array  of  N  +  1  point  sensors  with  uniform  separation  lA.  Let  the  total  length  of  the  array  be 
La  -  NIa .  We  fix  time  at  some  instant  and  decompose  the  spatially  fluctuating  field  80  (x)  into  a  spa¬ 
tial  spectrum.  In  accordance  with  the  Nyquist  limit,  all  spatial  periodicities  less  than  21  A  will  not  be 
effectively  measured.  The  distance  2 lA  is  thus  the  minimum  distance  of  response  to  the  spatially 
fluctuating  field.  The  upper  cutoff  wavenumber  is,  therefore, 

“max 355  -7“  (units:  rad/m). 

uA 

Of  the  remaining  wavenumber  band,  0  to  amax,  a  second  portion  is  cut  off  by  the  length  of  the  total 
field  sensed  LA  (i.e.,  track  length),  and 

“min  ^  ~r~  (units:  rad/m). 

Between  the  two  cutoffs  one  then  finds  the  effective  wavenumber  band  of  the  spectrum  to  be 

A  11  ( _ 


At*  “min  “max 


(units:  rad/m). 


If  the  array  is  continuous  rather  than  a  collection  of  point  sensors,  the  band  is 

A  __  1  ( _ _ 1/ _ \ 


Aa  =  oo  — 


(units:  rad/m). 


Instrument  windows  have  their  own  spatial  and  temporal  spectra.  Let  w(x,t)  be  the  window 
modeled  as  the  impulse  response  of  a  filter.  The  filter  transfer  function  is  then 


W(ay  w)  —  J*  w(x,  t)e 


iwt—iax  dtda 

(2t t)2 


The  variance  of  a  random  scalar  field  (say  temperature)  is  then  given  by 


<802>  —  J  J  1^80 (“» <o)  |2  |  IV (a,  (o)  |2  da  dw. 


/ 
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Here  the  spatial  and  temporal  limits  imposed  by  the  window  permit  the  integration  limits  to  be  ±oo. 
The  window  effectively  cuts  off  the  contributions  of  the  a—oj  spectrum  of  the  fluctuating  field  when 
these  extend  beyond  its  own  limits. 


Closed  Form  Solutions  for  Mean-Square  Temperature  Fluctuation 


For  the  remainder  of  this  report  we  consider  in  detail  solutions  of  the  basic  equations  for  band- 
limited  mean-square  temperature  fluctuations  in  the  ocean.  It  was  noted  earlier  that  for  a  Gaussian 
autocorrelation  of  the  temperature  field  O(xj)  the  band-limited  variance  is 


<802( Aa,  A oj)>  « 


// 

AttAtu 


<802>  lcUTcfi  -- 
2n  e 


dado) 


(units: 


The  integration  is  elementary  if  the  assumption  is  made  that  there  is  no  dispersion: 


where 


<802(Aa,  Acu)>  = 


<802> 


erfl 


Jcfi 


V2 


erf] 


V2 


x 


erfl 


V2 


a)minTc« 

V2 


erf(z)  =  -j=r 

V7T 


This  is  a  closed-form  solution.  Here,  as  before,  the  values  of  amax,  amin,  a)max,  and  o.»min  must  be 
determined  by  considerations  of  the  instrumental  and  experimental  cutoffs  to  the  a-co  spectrum 
described  earlier  in  this  report. 


We  next  consider  a  towed  sensor.  Let  the  speed  of  tow  be  v(m/s).  If  this  is  fast  enough  (say 
greater  than  a  few  knots)  the  ocean  temperature  field  is  essentially  "frozen."  Then  we  can  convert  the 
o>  spectrum  to  an  a  spectrum  by  writing  a  —  co/v.  For  a  Gaussian  autocorrelation  of  the  temperature 
field  we  find  the  closed  form  solution 


<802(Aco/v)> 


<6  92> 
2 


yj~2  V 


—  erf 


yflv 


These  band-limited  variances  of  the  temperature  field  will  be  used  later  in  making  numerical  estimates. 
This  completes  the  first  major  portion  of  our  investigation. 


We  now  take  up  the  second  major  problem,  which  is  to  make  numerical  estimates  of  the  effects  of 
temperature  fluctuations  on  the  phase  stability  of  optical-fiber  hydrophones. 


TEMPERATURE  FIELDS  IN  THE  OCEAN 


The  temperature  field  0(r,/)  in  the  ocean  is  a  complex  phenomenon.  It  consists  of  several  com¬ 
ponent  fields  distinguished  among  themselves  by  differing  length  and  time  scales.  First,  there  is  a  long- 
duration  mean-temperature  variation  with  range  and  depth,  labeled  the  temperature  profile.  Secondly, 
there  is  a  faster  space-time  fluctuation  of  temperature  at  each  range  and  depth  about  the  mean  profile. 
Thirdly,  there  is  a  zi\\\  shorter  scale  space-  and  time-varying  microstructure  of  temperature  which  also 
varies  with  range  and  depth.  The  causes  of  these  component  phenomena  are  thought  to  be  seasonal 
solar  heating  which  forms  the  profiles,  internal  waves  generated  by  winds,  convection  currents,  etc. 
which  distort  the  profiles  in  space  and  time;  and  velocity  turbulence  and  thermal  diffusion,  which  are 
responsible  for  creating  the  temperature  microstructure. 
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We  concentrate  here  first  on  the  fluctuations  in  temperature  profile  caused  by  internal  waves  and 
second  on  the  temperature  microstructure.  Our  objective  is  to  calculate  band-limited  variances  of  tem¬ 
perature  fluctuation  in  each  category  and  from  these  to  calculate  the  mean-square  fluctuation  in  phase 
of  a  laser  light  beam  being  transmitted  through  a  fiber  submerged  in  the  random  temperature  field. 


We  begin  with  the  effect  of  internal  waves  on  the  temperature  field.  Before  entering  the  analysis, 
we  define  the  concept  of  one-dimensional  spectra,  which  are  the  major  tools  for  calculation. 


One-Dimensional  Spectra 


All  spectra  used  in  this  report  are  one  dimensional.  They  are  derived  from  three-  or  four- 
dimensional  spectra  by  integration  in  the  usual  way.  Let 

<80(xj)  89(x  +  t  t  +  t)>  =  <892>  p(£.  r> 

=  f~  ,OJT  dadw, 

—  oo 


where  p  is  the  nondimensional  correlation  function  and  8S  is  the  spectral  density.  Then  the  one¬ 
dimensional  spectra  are  generated  as  follows: 


<892>  p(a j.  0.  0,  0)  —  J  ^  e  *rt|  8SH(a,  a»)  dadio 

=  J  ^  J J'J*  SSH(a,  (jj)  da2da^dw 

— oo 

=  f°°  p'*1"1  F'(a, )  da,. 

*/  -  oo 


(1) 


da  1 


where  Fl(a\)  is  the  one-dimensional  spectral  density  along  coordinate  and 

<892>  =  p(0,  0,  0,0))=  f  e~><OT  F*(aj)  dio, 

•J  —oo 


oo 

F*((o)  =  Iff  e,a *  85 (a,  (o )  da, 

—  oo 


(2) 


where  F*( a>)  is  the  one-dimensional  spectral  density  along  the  coordinate  of  time.  If  the  field  80  is  iso¬ 
tropic  in  space,  then 

8  S  (a ,  co )  =  8S(a,oj),  a2  —  a  *  a. 


If  the  field  is  not  isotropic,  it  is  sometimes  convenient  to  define  a  three-dimensional  spectrum  of  a  sin¬ 
gle  variable  by  integration  over  a  spherical  shell  of  radius  a: 

E(a,io)  =  si  8S  (a,a>)  dA. 

a  a=a‘ 


Then 


J»  oo 

E(a)  da. 

o 


Garrett  and  IMunk  Spectrum  for  Internal  Waves 

Garrett  and  Munk  proposed  a  model,  called  GM72  [3],  later  improved  to  GM75  [6],  of  a  power 
spectrum  of  vertical  displacement  of  the  ocean  medium  at  various  depths  caused  by  internal  waves. 
The  model  is  given  by  their  Eqs.  (6.20)  and  (6.22): 

n  F^iio)  —  Ewi  J-r  ~~ (or  -  u)2)^ 

l  7 7'  M  it) 


(GM6.20) 
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h  Ft(a,)  =  y 


1  r  I  c,  1 

— r  Eu>,  N  — , 

ttj  M3  a,2 


«  «  n. 


(GM6.22) 


The  symbols  used  by  them  are  very  special.  They  are  listed  in  Table  1  together  with  their  units. 


Table  1  —Garrett  and  Munk  Symbols 


Symbol* 

Meaning 

SI 

Units 

GM 

Physical  Units 

h 

Vaisala  frequency 

s^ 1 
nr  s 

cph 

Fi*») 

E 

G>( 

vertical  displacement  oj  spectrum 
nondimensional  energy 
nondimensional  inertial  frequency 

m2/cph 

M 

(U 

scale  factor  for  spatial  coordinates 
nondimensional  frequency 

rrf1 

cpkm 

vertical  displacement  a  spectrum 

m3 

m2  •  cycle/m 

N 

V 

scale  factor  for  time  coordinate 
nondimensional  wavenumber 

s'1 

cph 

d 

wavenumber 

m"1 

cpm 

inertial  frequency 

S'  i 

m3  •  s  2 

cph 

E 

energy  in  units  of  velocity2  x  distance 

p 

mass  density 

N  ■  s2  *  rtf 4 

pe 

energy  per  unit  area 

N  •  m/  m2 

J/cm2 

Ely) 

energy  in  units  of  velocity2 

m2  *  s  2 

N  *  m/  m3 

J/cm3 

pE(y) 

energy  per  unit  volume 

x.y 

range,  depth 

— buoyancy  depth  scale 

2t tM 

m 

b 

m 

m/rad 

•Also,  a  =  a/M*  w  *  u/N,  n  -  w/A,  x  -  2nMx%  t  -  2nNt,  and  an  exponential  model  of  h  is 
given  by  h  *=  where  b  ~  1.3  km  and  y  >  200m. 


By  curve-fitting  to  ocean  data  GM72  assigned  the  following  values: 

M  =  0.1 22  cycles/km -1.22  x  l(T6  cycles/cm  =  1.22  x  10  4  cycles/m. 

AT1  =  8.196  x  105  cm/cycle  =  8.196  x  101  m/cycle, 

N  =  3  cycles/h  =  8.333  x  10"'4  cycles/s. 
lCr}N2  =  (8.196  X  105)’  (8.333  x  10  4)2  =  3.824  x  10"  cnrVs2. 

E  =  3.82  x  lO6  cnrVs2 . 
f,E  =  0.382  J/cm2  =  3.82  x  106  erg/ cm2. 
p  =  1  g/cm1  =  1  dyne  ■  s2/cm4  (1  kN  •  s2/m4). 

E  =  Ehr  =  3.82  x  IQ1'  x  2n  =  2ir  x  10  5 
XT2N2b  \fyN2  3.82  x10" 

J,,  =  0.04 cph  =  1.111  x  10 cycle/s  (at  30°lat.), 
j,  =  Jk  =  20  modes,  and 
b  =  1.3  km  -  1.3  x  10' m  =  1.3  x  10s  cm. 

In  order  to  obtain  a  quantitative  feeling  of  these  entities  we  cite  a  few  examples. 


Example  /:  In  Fig.  2  (GM72  [3]  Fig.  4)  choose  a,  =  1  cpkm  (cycle  per  kilometer),  then 


1  cpkm 
0.122  cpkm 


8.196 
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>  s.  i* v — 


and 


m«,)= 


6.64  x  lO’mVh 
2.38  cph 


=  2.79  x  103 


nr 


cycle 


or 


m 


cpm 


This  also  checks  with  the  value  found  in  Fig.  2. 


Example  3:  According  to  GM72  Eq.  (6.20), 

»''/•„  (<*>)  =  —  £«(ArV3 

7 T 

Now  choose  o>  =  0.1  cph, 

0_l_C£h  =  3333 x  2 

N  3  cph 
and 


co2  +  <o 2 


a>2  —  cu 2 


units: 


m 


cycle2 


Thus, 


h 


03, 


O04^h  _  ,  333x  io-2 
N  3  cph 


/■„<«)-  — 

77 


x  (2ttx  10~5)(1.333x  1Q~2)  (8.I96x  IQ3)2 
(3.333  x  10“2)3 


(3.333  x  IQ  2)-  +  (1.333  x  IQ  2)2 
V (3.333  x  10"2)2  -  (1.333  x  1(T2)2 


=  4.08  x  104  m2/ cycle2. 


This  checks  with  the  value  shown  in  Fig.  2. 

The  proposed  GM72  model  also  gives  a  relation  between  vertical  displacement  spectra  and  vertical 
velocity  spectra: 


Thus,  if  oj2  »  (o 2,  then 


n  {  Fu( <d)  j  o>2  +  032 

— .--X  =  4 77*  — - r 

n  ((d)  o) *  o> / 

n~{  Fu (03)  — ■  h  Fr((o)  477 2 


or 


Fu(uj)  — • *  n2  F^(cd)  477 2 

One  can  then  find  /^(o;),  if  convenient,  then  obtain  F^ico). 
Ocean  Vertical  Displacement  Spectra 


units: 


mij_ 

h2  cph 


Power  spectra  of  internal  waves  are  classified  by  GM72  according  to  the  method  of  measurement: 
towed  spectra  (TS),  obtained  by  towing  temperature-sensitive  sensors  along  a  horizonal  line  in  the 
ocean;  moored  spectra  (MS),  obtained  by  fixing  the  sensors  in  a  location  in  the  ocean;  or  dropped  spec¬ 
tra  (DS),  obtained  by  dropping  weighted  sensors  in  the  ocean.  In  each  case  the  Garrett  and  Munk 
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models  of  spectra  in  the  three  spectral  coordinates  (horizontal  wavenumber  au  vertical  wavenumber 
a 2,  and  frequency  o>)  are  collapsed  by  integration  into  a  single  coordinate.  The  selected  single  coordi¬ 
nates  are  a\,  in  TS;  w,  in  MS;  and  a2 ,  in  DS.  Examples  of  spectra  used  in  this  report  are  found  in 
GM72  [3]  and  GM75  [6].  Figures  3a,  3b,  and  3c  reproduce  Figs.  2,  3,  and  5  of  GM75. 

In  each  of  these  experimentally  obtained  spectra  the  corresponding  models  of  GM72  and  GM75 
are  superimposed.  Of  particular  importance  are  the  upper  limits  of  measurement:  ai  =  3xl0~2 
cycle/m  in  TS,  <d  =  8  cycle/h  in  MS,  and  find n  =  30  cycle/m  in  DS;  and  the  lower  limits  of  measure¬ 
ment:  aj  =  3x  10~5  cycle/m  in  TS,  a>  ~  10“ 1  cycle/h  in  MS,  and  pnjn  =  3x  1CT2  cycle/m  in  DS. 
These  limits  impose  restrictions  on  the  possible  cutoffs  in  the  calculation  of  band-limited  variances 
undertaken  later  in  this  report. 

Figure  3a  (GM75  Fig.  2)  also  shows  a  spectrum  of  moored  vertical  coherence  (MVC).  This 
report  does  not  discuss  explicitly  coherence  and  its  effect  on  noise  calculations.  We  intend  to  investi¬ 
gate  these  effects  in  subsequent  reports. 


Fig.  3a  —  Moored  spectrum  in  coordinate  iu  obtained 
by  J  L.  Cairns  |J.  Geophys.  Res.  80,  299-306  (1975)1 
(GM75  (61,  Fig.  2) 


Fig.  3b  —  Improved  normalized  towed  spec¬ 
trum  in  coordinate  a,  made  by  E.J.  Katz  IJ. 
Phys.  Oceanogr.  3,  448-457  (1973))  (GM75 
(6) ,  Fig.  3) 
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Fig.  3c  —  Normalized  dropped  spectrum  in  coordinate  {nj n)($  — 
( n0/n)a2  obtained  by  R.  Millard  [MODE  Hot  Line  News,  No.  18,  1 
(1972)]  (GM75  [6],  Fig.  5) 


Effects  of  Tow  Speed  on  the  Spectrum 

The  a|,  a2,  oj  spectrum  of  vertical  displacement  of  internal  waves  describes  both  spatial  and  tem¬ 
poral  aspects  of  the  vertical  displacement  as  natural  phenomena  of  the  ocean.  In  contrast,  the  speed  of 
tow  describes  how  fast  the  phenomena  are  sampled  in  time.  At " infinite”  speed  the  sampling  is  instan¬ 
taneous.  Hence  the  resulting  spectrum  is  that  of  a  “frozen  ocean."  It  is  a  wavenumber  spectrum.  At 
zero  speed  of  tow  the  sampling  in  spatial  terms  vanishes,  leaving  only  a  sampling  in  time.  The  data 
contribute  only  to  a  frequency  spectrum.  The  model  of  these  spectra  used  in  this  report  is  that  of  Gar¬ 
rett  and  Munk,  GM72.  This  model  specifies  the  zero-speed  spectrum  as  the  moored  spectrum  (MS)  and 
the  infinite-speed  spectrum  as  the  towed  spectrum  (TS).  At  intermediate  tow  speeds  the  sampling  time 
is  slow  enough  to  permit  the  ocean  to  vary  somewhat  in  time  during  sampling.  GM72  models  these 
intermediate  speeds  in  Fig.  9  of  their  article,  reproduced  here  as  Fig.  4.  At  low  speeds  the  spectra  are 
plotted  as  hF^icr)  vs  o-,  where  o*  =  |w  —  aS\  is  the  "frequency  of  encounter."  They  are  essentially 
modified  moon  d  spectra.  At  high  speeds  the  spectra  are  plotted  as  hSF^(a\)  vs  a/S.  They  are  essen¬ 
tially  wavenumber  spectra  which  merge  into  the  GM72  model  of  towed  spectra  when  the  speed  is 
infinite.  Figure  4  (GM72  Fig.9)  also  shows  that  intermediate  speeds  strongly  affect  only  the  upper  and 
lower  limits  of  the  spectra.  The  influence  of  tow  speed  on  the  middle  ranges  of  the  spectra  is  minimal. 
Thus,  for  ordinary  tow  speeds  (say  2  knots  and  above)  the  ocean  is  essentially  frozen  for  internal 
waves,  except  for  the  very  shortest  and  very  longest  wavelengths.  However,  actual  experimental  data 
(see  Fig.  2  [Fig.  4  of  GM72))  do  not  agree  with  the  model  at  the  upper  and  lower  extremities  of  the 
spectra.  These  discrepancies  have  been  removed  in  a  later  model,  GM75  [6].  The  calculations  of  this 
report  are  unaffected  by  them,  since  they  are  based  on  the  middle  ranges  of  the  spectra,  the  model  of 
which  is  not  changed  in  GM75. 


Calculation  of  the  Mean-Square  Fluctuation  in  Temperature  over  a  Towed  Fiber 

A  fiber  of  length  L  —  1000  m  is  towed  at  a  depth  of  v  —  100  m  at  a  speed  S  knots  such  that  it 
floats  horizontally.  We  wish  to  calculate  the  mean-square  fluctuation  in  temperature  <802>  over  the 
entire  length  for  all  temporal  frequencies  and  vertical  wavenumbers. 
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f-  J 

b 

^  103 


01  <>,  cph  1 


ioH 


O.OI  0.1  i  10  o.oi  0.1  1  10 

o\  cph  o7 S,  cpkm 

Fig.  4  —  Doppler  effects  on  the  current  spectrum  (left)  and  vertical  displacement  spectrum 
(below),  for  n  =  1  (3  cph),  E  —  2n  x  10“ 5  (0.38  joules/cm2),  jt  —  20.  For  low  speeds  the  spectra 
are  plotted  against,  frequency-of-encounter  a\  for  high  speeds,  against  ar/S.  The  curves  are  for 
relative  speeds  S/S0  *  0,  0.45,  0.90,  1.80,  <»(S  -  0,  10,  20,  40,  °°  cm/ sec)  with  0,  oo  correspond¬ 
ing  to  the  moored  and  towed  spectra  (6.20,  6.21)  respectively.  (GM72  [3],  Fig.  9) 

To  make  this  calculation  we  shall  use  the  convenient  formula  [3,  footnote  on  page  237]: 

<802(Aaj,  _y)>  =  -^7-  f  Fc(a\ty)da\=  —  <i2> 

dy  •'A®  1  dy 

Thus,  we  must  first  find  the  mean-square  vertical  displacement  <f2>  and  the  temperature  gradient 
dT/dy  at  the  depth  y  in  question.  To  do  this  we  require  ocean  data  and  a  suitable  model  of  the  power 
spectrum  Fc(a\,  y).  A  well-substantiated  model  has  been  proposed  in  several  articles  by  Garrett  and 
Munk  [3,6].  We  assume  in  what  follows  that  ^(<*1)  is  the  same  as  in  GM72  [3]  Eq.  (6.22).  This  is 
called  by  them  the  towed  spectrum ,  and  is  a  model  developed  from  the  data  of  many  towed  sensor 
experiments.  It  includes  the  effects  of  all  temporal  frequencies  and  all  vertical  wavenumbers  associated 
with  internal  waves.  A  more  elaborate  model  of  towed  spectra  is  found  in  Fig.  4  (Fig.  9  of  GM72).  It 
shows  that  the  parameter  S  (tow  speed)  affects  the  spectrum  appreciably  only  at  the  upper  and  lower 
ends  and  very  little  over  most  of  the  portion  in  between.  We  shall  use  here  only  the  simplest  GM72 
version  of  the  tow  spectrum,  which  has  the  form 


,  1  1  .  1 
(«»)  ~ - -  /,  — r 


E  cu./V 


o>.  <  co..  <  n. 
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The  definitions  and  values  assigned  to  the  symbols  are  taken  from  GM72  [3]  and  Bell  et  al.  (7j: 


h 

Ji 

E 

*>/ 


dT 

dy 

M 

N 

«i 


=  Vaisala  frequency— at  100  m,  h  =  5.0  cycle/h; 

=  number  of  vertical  modes— GM72  takes  j,  =  20; 

=  dimensionless  energy  constant— GM72  assigns  E  =  x  10~5; 
—  dimensionless  inertial  frequency— GM72  assigns 
004^  0_2 

N  3cph 

=  temperature  gradient  at  depth  y 


=  0.025  at  100-m  depth  (a  discussion  of  temperature 
gradient  is  undertaken  in  the  next  section  of  this  report); 
=  spatial  scale  factor— GM72  assigns  M  =  0.122  cycle/km; 
=  temporal  scale  factor— GM72  assigns  N  =  3  cycle/h; 

=  dimensionless  wavenumber— GM72  assigns 


M  0.122 


8.196a!  , 


where  a  i  is  the  wavenumber  in  units  of  cycle/km. 


We  calculate  the  constant  factors  of  Fc  in  units  of  meters  and  hours: 

i  1  •  J_  Ea>N  =  _LxixinJ_y  2tt  x  IQ- 5  x  1,333  x  IQ-2  x  3 
ft  2  J'  it*  M3  5.0  2  n3  '  (1.22  x  10~4)3 

=  8.925  x  104  m3/cycle  (=m2/cpm). 


The  integral  required  is 


a  ry  7-  J  Jl? 


da  i 


«lr 


<*ln 


Since  the  upper  portions  of  the  spectrum  contain  little  energy,  we  arbitrarily  take  The 

minimum  wavenumber  must  be  estimated.  If  we  assume  slow  tow,  and  the  effective  length  of  the  data 
sample  is  1000  m,  we  estimate 


<*ln 


__  1  cycle  _  1(v-3 


1000m 


lO-3  cpm. 


Thus, 


and 


/  =  (1.22  x  10-4)2  x  =  1.488  x  1 0'5  cpm 


10 


<£2>  =  f  Fr  (a,)dat  =  8.925  x  104  x  1.488  x  10“ 5  cpm 
J  *  cpm 


=  1.33  m2. 

This  means  that  over  a  length  of  1000  m  the  rms  vertical  displacement  at  a  depth  of  100  m  is  about  1 
m.  The  rms  temperature  fluctuation  over  the  length  of  the  fiber  is  thus 

<d92(Sa,y)>  -  (0.025)2  x  1.33  -  8.3  x  10”4  (°C)2. 

Note  that  this  result  includes  the  effect  of  all  frequencies  [6]. 
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Discussion :  Clearly  at  a  different  depth  both  the  Vaisala  frequency  and  the  temperature  gradient 
will  be  different  from  what  is  assumed  above.  Thus  the  calculation  will  be  strongly  depth  dependent. 
Also  it  must  be  recognized  that  GM72  does  not  give  a  good  approximation  if  the  tow  depth  is  less  than 
100  m.  In  these  shallow  depths  the  mixing  of  different  temperature  strata  partially  destroys  the  internal 
wave  structure. 


To  check  our  result,  we  asked:  What  will  be  the  value  of  <£2>  if  the  entire  spectrum  of  energy 
is  considered?  To  answer  this  we  made  the  following  calculation. 


Let  us  consider  the  case  of  towed  spectra  given  by  Fig.  5  (Fig.  1  of  Bell  et  al.  [7]).  We  can 
approximate  the  function  Fc  (c*i)  (units:  m2/cpm)  by  the  formula 

-  t-2 


Fc(a{)=  105 


<*i 


10“ 


105  x  10"8  1(T3 


a  i2 


«,2 


1 


1 


in  which  the  units  of  a,  are  cpm.  To  check,  take  a|  =  10  3  cpm;  then  Ft  ~  103,  which  agrees  with 
Fig.  5  (Fig.  1  of  [7]).  Now 

dot  ] 

Choose  at  \ max  — *  00 ;  then 

/  =  10“3 


<C2>  =  J  F(  («,)</«,  -  10“3/  -  1(T3 


«ln 


a  lmin 


Assume  a  |min  =  10  4  cpm;  then 


1 


I  min 


<£2>  =  10m2  and  V <£2>  —  3m. 
For  the  range  shown  in  Fig.  5  (Fig.  1  of  [7]), 


<£2>  =  10-3 


i 


I 


10-4  10- 


=  10  m2. 


Thus  at  a  depth  of  100  m,  a  towed  sensor  will  show  a  root-mean-square  value  of  vertical  displacement 
of  about  3  m  when  all  wavenumbers  from  10-4  to  00  cpm  are  considered.  In  the  earlier  calculation  we 
considered  the  spectrum  contributions  at  the  lower  end  to  stop  at  10~3  cpm,  which  gave  the  result 
<£?>  =  1.33  m2.  The  corresponding  temperature  fluctuation  for  the  entire  spectrum  is  <502>  =  6.2 
x  10' 3  (°C)2.  Note  again  that  the  spectrum  (aj)  includes  the  effect  of  all  frequencies,  as  required 
both  by  definition  and  construction. 


Calculation  of  Phase  Fluctuations  Caused  by  Temperature  Fluctuations 


We  noted  in  the  section  "Environmental  Noise"  that  the  phase  fluctuation  80  due  to  temperature 
fluctuation  50  (x)  with  correlation  length  Lc  in  a  fiber  of  length  £,  of  index  of  refraction  w,  transmitting 
a  laser  light  of  wavelength  X0>  is 


•\/  <50/  >  =  C^£*\/  2  LLC  <802> . 

Assuming  that  ILL,  =  0(1)  m2  and  using  the  explicit  form  of  C9L ,  we  write 


V  <M>[>  =  k0n 


1  8n  1  9L 
n  d6  L  d8 


V<8»2>. 


Let  us  choose  a  glass  fiber  and  make  the  additional  specifications 

n  =  1.46, 

X0-  5.6  x  103  A. 
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Fig.  5  —  Normalized  horizontal  (A)  and  vertical  (B)  one-dimensional  spectra 
of  internal  wave  amplitude  (Bell  et  al.  (7],  Fig.  1) 


and 


Then 


and 


Hence, 


1  4 

L  89 


x 


icr7 


units: 


_1_ 

•c 


|f-  11.9  x  10-6 
oO 


units: 


V 


*n  = 


2tt 


2ir 


5.6  x  10  7 m 


|7 


—  =  1.122  x  10' m 


k0n\ 


T|f  +  1  =  1.122  x  107 

H.9  x  10-6  ,  4  x  ,o~7l 

n  89  L  89  J 

l 

1.46 

J 

=  1.40  x  102 

rad 

units: 

‘C  •  m 

x  1.46 


V<6^7>  =  1.40  x  102  V<802>. 


We  continue  the  calculation  of  -J <S<^2>  using  the  GM72  spectrum  for  various  fiber  lengths  L.  The 
results  are  shown  in  Table  2.  Actual  calculations  are  shown  in  Appendix  F. 
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Table  2— Values  of  V  <8 <£/> 
for  Various  Fiber  Lengths  L 


L 

V<8  <f>2> 

(m) 

(rad) 

103 

4.04 

102 

1.28 

10 

0.404* 

1 

0.128* 

‘Extrapolated  values 
(uncertain) 


Temperature  Gradient 

The  temperature  gradient  dT/dz  in  the  ocean  is  a  function  of  depth,  season,  salinity,  etc.  In  any 
particular  calculation  of  temperature  fluctuation  it  is  essential  to  measure  this  gradient  at  the  locale  in 
question.  However,  it  is  equally  useful  to  have  a  model  of  this  parameter. 


A  convenient  starting  point  in  the  development  of  this  model  is  the  equation  of  sound  speed  gra¬ 
dient  proposed  by  Munk  [8]: 


where 


9*-C  dr.„dS 


+  7a> 


a  =  3.16  x  HT3^)-1; 

0  =  0.96  x  10"3(%)-1; 
y A  =  1.14  x  1(T2  krrr1; 
d:C  =  dC/dz\ 

C{z)  is  the  speed  of  sound  as  a  function  of  depth; 

T  is  the  potential  temperature,  which  is  the  observed  temperature  minus  the  adiabatic  gradient; 
S  is  the  salinity;  and 

yA  is  the  fractional  velocity  gradient  in  an  adiabatic  ocean. 


We  first  eliminate  salinity  effects  by  use  of  the  Turner  number ,  Tu: 


Thus, 

and 

where 


dS  _  a  dT 
dz  b  dz 


Tu. 


9.-C  dr  ,  8a  BT  -  . 

~TT  =  a  -r—  +  t~  -r—  Tu  + 


C 


dz 


b  dz 


7  a 


dr  _  £%c  -  -ya 

dz  a(l  +  cTu) 


C  s= 


a§_ 

ab 
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The  fractional  sound  speed  gradient  is  related  to  the  Vaisala  frequency  N(z).  In  Munk’s  model, 


C~ldzC 


JL 

g 


N2(z)  +  yA, 


and 


Thus, 


/x  =  —  s{Tu), 
a 

(rr  v  1  +  cTu 

s(Tu)  =  - - — , 

1  -  Tu 

—  =  24.3, 
a 

c  =  =  0.049. 

ab 


1L 

N2(z) 

dr  _ 

g 

a(l  +  cTu)  ..  1 

NHz) 

dz 

a  < 

(1  +  cTu)  tf£(l  —  Tu)  a 

(1  +  cTu) 

NHz) 

ag(l  -  Tu)  ' 

A  convenient  model  of  the  Vaisala  frequency  below  the  thermocline  is  the  exponential  one, 

N  =  N0e~z/B, 

N0=  5.24  x  10~3—  (=  3  cph), 
s 


and 


B  —  1.3  km. 


In  the  Atlantic  near  Bermuda  Tu  is  about  +0.8.  Thus,  with  depth  in  meters, 


ar 

dz 


(5.24  X  lQ~3)2e-2z/1.3».103 

1.3  X  10-4  X  9.8(1  -  0.8) 


(units:  lC/m). 


Hence,  at  z  —  100  m, 


2xl02 

or  _  (5.24  X  10~3);  e  1 3xl°3 

dz  1.3  x  lr4  x  9.8(0.2) 


=  -9.3  x  10~2  “C/m. 


In  the  North  Pacific  Tu  =  -0.3.  Hence  at  100-m  depth  —  — 1.4  x  10-2  °C/m.  These  are  average 

oz 

values  over  whole  ocean  basins. 


Vertical  Temperature  Fluctuations  at  Fixed  Depths 

Bell  et  al.  [7]  lowered  70  thermistors  in  a  vertical  chain  at  1.4-m  separation  to  a  mean  depth  of 
100  m  and  towed  them  at  slow  speed.  From  the  measured  data  they  constructed  a  one-dimensional 
(vertical)  spectrum  of  internal  wave  amplitude.  This  is  shown  in  Fig.  5b  (their  Fig.  IB).  The  abscissa 
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is  (hjh)  cpm,  and  the  ordinate  (/?/ w0)2 in  units  of  m2/cpm.  It  is  important  to  note  that  the 
ordinate  includes  a  (natural)  integration  over  all  horizontal  wavenumbers  and  all  frequencies  which  lie 
within  the  range  of  the  length  of  the  tow  track  and  the  duration  of  the  experiment. 

For  purposes  of  making  estimates  we  approximate  the  functional  dependence  of  Bell  et  al.  (Fig. 
5b)  by  a  simple  power  dependence. 


y  =  ( n/fto )  F^ip) 


n°  « 

*  =  —  0* 
n 


3.184  x  10~2 


Now  the  integral  over  the  entire  spectrum  from  x  =  10  to  x  -  10  3  is 


L ydx  -  L 


3.184 x  10'2  J  3.184  x  10“2 
- - dx  — - — - x  11 


-3,184  x  1Q~2  1 

1.2  (10)1-2 


do-3)1 


=  106. 

However,  the  range  of  wavenumbers  resolved  by  the  experiment  of  Bell  et  al.  was  smaller.  At  the 
upper  limit  the  array  sampled  the  vertical  temperature  field  at  an  effective  spacing  of  1.27  m  between 
thermistors.  The  nominal  limiting  wavenumber  was  1 . 2 7“ 1  =  0.787,  and  the  Nyquest  limit  was  half  of 
this,  or  3.9  x  10“’  cpm  (say  4  x  10'1  cpm).  The  lower  limit  in  the  vertical  spectrum  was  obtained  by 
using  16  decrements  from  the  upper  limit  of  size  4  x  1 0”  V 1 6  =  2.5  x  10~2  cpm,  so  that  the  lowest 
wavenumber  resolved  was  2.5  x  10" 2  cpm,  corresponding  to  a  physical  length  of  2/ (2.5  x  10"2),  or  80 
m.  The  range  of  vertical  wavenumbers  was  then 

2.5  x  10“ 2  <  p  <  4  x  10“l  cpm. 

Now  at  a  depth  of  100  m  the  Vaisala  frequency  was  measured  at  —  5.0  cph.  Hence  n/n0  ~  5/3  = 
1.67.  Thus  the  range  of  parameter  x  =  ( n/h0)p  was 

4.2  x  10~2  <  x  <  6.7  x  10“!  cpm. 

The  mean-square  displacement  in  the  range  was  then 

(  .  U  f  .  i  _ _  .  _  r  i 


r  -  )2  . 

<£2>  =  J  —  F£(/3)  d  —  i 3  = 

l  no)  «o 


-3.184  x  ir2 

1.2 

1.148=  1.2  m2. 


(6.7  x  1 0“ 1 ) 1 2  (4.2  xlO-2)1 2 


Thus  the  mean-square  vertical  displacement  over  a  vertical  length  of  nominally  100  m  centered  at  a 
depth  of  100  m  is  —  1.2  m2.  The  mean-square  fluctuation  in  temperature  is  then  calculated  from  the 
measured  value  of  the  temperature  gradient  dO/dy  at  that  depth: 


*  ✓ 


~~  »- 
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<802(A/3,  y)>  = 


(W 

dy 


<C> 


=  (0.025)2  x  1.2  =  7.5  x  Kr4m2. 


This  does  not  differ  materially  from  the  mean-square  temperature  fluctuations  calculated  from  the  hdr- 
izontal  spectra  noted  above.  Using  the  same  numerical  values  of  k  and  n  as  before,  we  obtain 

V<8^1>  =  1.40  x  102  x  VT5  x  10'4  -  3.83  rad. 


Calculation  of  Temperature  Fluctuations  for  Moored  Spectra 

Assume  that  a  single  point-source  sensor  is  moored  at  a  depth  of  100  m,  where  it  detects  changes 
of  temperature  with  time.  A  composite  model  of  the  spectrum  of  vertical  displacement  at  this  depth  is 
given  by  GM72  [3].  We  use  this  moored  spectrum  (MS)  to  calculate  the  mean-square  fluctuation  in 
temperature.  Thus,  the  mean-squared  vertical  displacement  is 

<d2)1/2 


i  C  *  ,  ,  111  r  (<°2 

<£2>  -  I  Fr  (oj)  -ry  I  - 

J  c  h  2  tt3  M2  J 


doj . 


The  definitions  and  values  assigned  to  the  symbols  are  taken  from  GM72  [3]  and  Bell  et  al.  [7]: 


h 

E 


M 

N 


=*  Vaisala  frequency— at  100  m  depth,  n  =  5.0  cph; 

=  number  of  vertical  modes— GM72  takes  j)  *  20; 

—  dimensionless  energy— GM72  assigns  E  =  2n  x  10~5; 

=  dimensionless  inertial  frequency— GM72  assigns 

M4ah_U33x  10-2 
N  3cph 

=  spatial  scale  factor— GM72  assigns  M  ~  1.22  x  10-4;  and 
=  temporal  scale  factor— GH72  assigns  N  —  3  cph. 


With  these  assignments, 


j  r  (cd2  —  w  2)l/r2 

<C2>  =  const  x  J  - - - - 


dco. 


Now, 


Next, 


n2i r3 


Eoj,  N 


Ixixix  —  -  10~5  x  1-333  x_10:2  x  M  =  05444m  2. 
5.0  2  7t3  (1.22  x  IQ"4)2  1 


f 


(fa)2  -  CO,2) 1/2 


d  CD 


f 


VRdx 


where 


R  =  a  +  bx  +  cx2  and 
A  =  4  ac  —  b2 . 


In  this  case  a  =  to,2,  b  =  0,  c  =  l,  A  =  -4<o,2. 


f  JRdx 

J  x3 


■JR  ,  I  f  dx 
lx2  2  J  xsfR 


Now  from  Gradshteyn  and  Ryzhik  [9.  p.  84], 


\Tr  1 

lx2  2 


arcsin 


2  (—to,2) 
Xyf—  4(— fa)2) 
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or 


c  (c2-*,,2)1'2 

•\f  OJ  2  ~  (O  } 

,  I 

oj, 

3 

*33 

3 

<1 

II 

2a,2 

T  — —  arcsin 
2(0 , 

(O 

V^max  ~  <*>,  1 

- — — -  +  - —  arcsin 

2 (o*  2a>, 


in  “W- 


2w 


2 


,  1  .  I  <*>, 

-4  - —  arcsin  | — 

2(o, 


^min  J 

Assume  that  the  entire  GM72  spectrum  is  under  consideration,  i.e.,  o>max  =  /)0,  wmin  =  d>;.  Now, 

o»n 


"'max  , 
w  max  “  —  1 

"0 


and 


Hence, 


‘"min  =  ~~  —  0P4C.Ph  =  1.333  X  IQ’2. 


«0 


3cph 


^  Su> 


Vl2-  (1.333  X  10“2)2 


1 


2  x  l2  2  x  1.333  x  I0“2 

V (1 .333  X  10"2)2  -  (1.333  x  10”2)2 


arcsin 


1.333  x  10'2 


2  x  (1.333  x  10“2)2 


+ 


1 


2  x  1.333  x  10 


r2 


1.0 

arcsin  (-1) 


=  +  1  -  [-0-  5.8919  X  10‘1  =  5.89  xio1. 

Thus,  the  mean-squared  displacement  for  all  frequencies  is 


and 


<{2>  =  0.5444  x  5.89  x  101  -  3.21  x  10 
>/<?>  =  5.7  m. 


m 


If  the  observation  time  is  restricted  to  be  less  than  w,  1  (*=  25  h/cycie),  then  the  lower  limit  of 
integration  is  increased.  Take,  for  example,  a  time  of  1  h.  Then  o>min  —  1  cph,  and 


;  =  — -min.  =  1  cph  =  3  33  x  10“ 1 
m,n  ho  3  cph  j  lU  • 


Hence, 


^ i\w 


-1 


arcsin  (— 1)  =  2.36, 


2  x  3.33  x  10"1 

so  that 

<{2>  -  0.544  x  2.36  =  1.28  m2. 

Finally,  if  the  observation  time  is  10  min,  we  take  wmin  ~  6  cph,  so  that 

_  6~ 
wmm  ^ 

I**,  =  Yx~2  arcsin  =  °-39, 
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and 


<{2>  =  0.544  x  0.39=  0.21  m2. 

The  mean-square  temperature  fluctuation  of  a  moored  sensor  at  a  depth  of  100  m  in  the  three  cases  is 
then  calculated  by  use  of  the  Bell  et  al.  temperature  gradient  of  0.025°C  [7]: 

For  a>min  =  a >, 

<S6>2>  =  32.1  x  (0.025)2  =  0.02  (°C)2 
and 

V<86»2>  =  0.14  °C. 

For  <umin  =  25(o, 


<802>  =  1.28  x  (0.025)2  =  8.0  x  10~4  (°C)2 


and 


For  wmin  =  1 50co  j 


and 


V<86>2>  =  2.8  x  10“ 2  V. 

<802>  =  0.21  x  (0.025)2  =  1.3  x  10~4  (°C)2 

V<802>  =  1.14  xl0“2  °C. 


To  find  the  corresponding  rms  phase  fluctuation  we  must  be  given  the  correlation  time  Tc,  the 
integration  time  T,  and  the  phase  conversion  factor  CgT.  In  the  absence  of  data,  we  assume 

2 TTC  =  0(1)  s2 


and 


IO,rl=  1  C*J=  1.40  x  102. 

Then  the  three  cases  become: 

-  W,  =  0.04  Cph, 

V<8#2>  =  1.4  x  102  x  0.14  =  19.6  rad; 

wmin  25(0,  “  1  cph. 

V<8#2>  =  1.4  x  102  x  2.8  x  10-2  =  3.92  rad; 

and 

w min  ~  150(i)/  =  6  Cph, 

<802>  =  1.4  x  102  x  1.14  x  10~2  =  1.60  rad. 

Thus,  if  we  allow  the  fiber  to  report  all  temperature  fluctuations  for  25  h  (1  cycle  ataj,  =  0,04  cph),  the 
rms  phase  fluctuation  is  very  large  (19.6  rad).  If,  on  the  other  hand,  the  cutoff  is  6  cycle/h  and  0  =  10 
min  observation,  then  the  rms  phase  becomes  1.6  rad. 

With  these  results  the  calculation  of  mean-square  vertical  displacement  of  internal  waves  from 
towed  spectra  (in  coordinates  a\  and  a2)  and  moored  spectra  (in  coordinate  &>)  has  been  completed  for 
particular  selections  of  depth,  range,  and  size  of  the  fiber-optic  sensor.  Corresponding  to  these  space- 
time  scales,  we  have  calculated  mean-square  temperature  fluctuations.  An  additional  calculation  can  be 
made  when  the  fiber  has  finite  mass,  and  hence  when  there  is  a  lag  before  steady  state  is  reached.  Let 
hit)  represent  the  difference  in  temperature  between  the  surface  of  the  fiber  and  the  center;  then 
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hit)  =  8(1  -  0)  -  j  e-'/r,  0<  /  <  T  (units:  s“‘) 

=»  0  outside  the  interval  0,  T, 


in  which  T  is  the  time  constant  of  the  equivalent  filter.  The  transfer  function  of  this  filter  is 


M(  a>)  = 


ijo  T 

1  +  iojT' 


whose  squared  modulus  is 


Thus  the  required  integral  is 


I//MI2 


(a >n2 

H(wT)2 


(to2  —  to2)1^2 


(toD2  ,  r  -Jr  dx 
[1  +  (qj  T)2]  J  Xl\  +0X2)' 


This  integral  can  be  evaluated  by  numerical  integration.  However,  for  values  of  T  such  that 
co2  T2  »  1  the  calculations  made  above  are  seen  to  be  valid.  Figure  6,  a  plot  of  I  // (cl» )  | 2 ,  shows  the 
effective  window. 


We  next  turn  to  vertical  displacements  (hence  temperature  fluctuations)  of  much  smaller  space- 
time  scales,  namely  the  microtemperature  field  of  the  ocean. 


Microtemperature  Structure  of  the  Ocean 

The  microtemperature  structure  of  the  ocean  has  been  measured  and  analyzed  at  an  increasing 
pace  in  the  last  15  years,  and  a  wealth  of  data  is  now  readily  accessible.  We  select  from  the  many  arti¬ 
cles  and  reports  a  few  examples  which  are  chosen  for  their  immediate  usefulness  in  making  numerical 
estimates  of  the  effects  of  microstructure  on  the  phase  stability  of  laser  light  in  a  fiber-optic  hydro¬ 
phone. 

We  begin  with  Gregg  and  Cox  1 1 01 .  They  measured  vertical  gradients  of  temperature  in  the 
ocean  on  a  millimeter  scale  in  the  depth  range  from  100  m  to  500  m.  The  site  of  measurement  was  the 
San  Diego  Trough.  The  method  used  consisted  of  a  free-fall  probe  with  no  connections  to  the  surface. 
Recovery  of  the  measurement  package  allowed  data  processing  of  the  tape-recorded  temperature 
records;  these  records  were  continuous.  The  analysis,  however,  was  based  on  a  sampling  procedure  in 
which  one  data  point  was  taken  every  1.5  mm.  The  mean  temperature  gradient  over  the  entire  100-  to 
500-m  record  was  +9.4  x  10~SoC/cm  (a  positive  gradient  means  temperature  increasing  upwards). 
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Chosen  for  analysis  were  two  sections  of  record  —7m  long  at  302-m  and  320-m  depth.  They  are 
reproduced  here  as  Figs.  7  and  8  from  Figs.  8  and  9  of  Ref.  10. 

We  discuss  Fig.  7  first  and  note: 

•  The  mean  vertical  temperature  gradient  over  the  record  was  -3.4  x  10“ 5  'C/cm. 

•  The  rms  gradient  (signs  omitted)  was  2.4  x  10~4  C/cm. 

•  A  net  inversion  of  temperature  of  magnitude  —  2  x  10-2C  occurred  over  a  vertical  thickness 
of  0.5  m  of  medium. 

•  Although  dO/dz  exhibited  random  structure,  (Hz)  was  relatively  smooth. 

Figure  8  exhibits  the  same  behavior  at  320  m  depth: 

•  The  mean  gradient  was  -2.2  x  10  s  °C/cm. 

•  The  rms  gradient  was  7.2  x  10~4  °C/cm. 

•  An  inversion  of  2.7  x  10'2  °C  occurred  in  about  2  m  of  vertical  thickness. 


302  metcfs 
SAN  DIEGO  TROUGH 
20  MAY  1970 


Z  -  1.0  metm 
(UP) 


TEMPERATURE 
GRADIENT 
- > 

~7’-  ♦3*l0*3oC/cm 


INTEGRATION 
OF  TEMPERATURE 
GRADIENT 

Hx10'2’C 


INTEGRATION 
OF  TEMPERATURE 
GRADIENT 


|^A«=  ♦  4  *10’*  *C 


Fig  7  —  Vertical  component  of  temperature  gradient  and 
its  integral  at  302  m  (10,  Eig.  8) 


Kig  8  —  Vertical  component  of  temperature  gradient  and 
its  integral  at  320  m  1 10,  Fig  9] 


Stommel  and  Fedorov  [11]  made  continuous  vertical  temperature  records  from  0  to  500-m  depth 
in  the  Phillipine  Sea  off  Mindanao  and  off  Timor  over  a  10-mile-square  range.  Figures  9a  and  9b  (Figs. 
3  and  4  of  Ref.  11)  show  the  salinity,  temperature,  and  depth  traces  of  Station  13  off  Timor  and  Station 
21  off  Mindanao.  Examination  of  these  traces  led  them  to  believe  that  sharp  vertical  gradients  of  tem¬ 
perature  about  2  to  40  m  in  thickness  existed  in  lamina  (or  stratification)  which  extended  horizontally 
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Fig.  9  —  Photographic  copy  of  actual  STD  trace  of  STD  station  13  off  Timor  [11,  Fig.  3)  (a)  and  STD  station  21  off  Mindanao 
(II,  Fig.  4]  (b).  The  depth  scale  is  500  meters,  each  small  division  being  5  meters.  The  salinity  curve  is  labeled  S,  the  scale  of 
salinity  is  labeled  at  two  points  34  and  35,  thus  each  smallest  division  in  the  abscissa  is  0.05 */•«.  The  temperature  curve  is  labeled 
T  (the  labels  A.  C,  and  P  refer  to  particular  points  of  interest  on  the  temperature  curve.  At  the  top  two  temperature  points  on 
the  abscissa  are  labeled  25°  and  26°C,  thus  each  smallest  division  is  0.05°C.  The  temperature  curve  has  origin  reset  at  successive 
depths  where  it  tends  to  go  off  scale.  This  is  done  in  steps  of  4 °C  as  indicated  by  the  numerals  2V,  17,  13,  9,  thus  point  A  is  at 
21.32°C.  The  pens  cannot  both  transverse  the  same  path  mechanically  so  the  temperature  curve  is  offset  downward  by  5  meters 
and  hence  always  reads  5  meters  too  deep.  On  tracings  this  has  been  allowed  for  and  corrected.  On  this  photograph  of  the  origi¬ 
nal  it  must  remain  uncorrected.  The  salinity  scale  is  not  offset  in  depth.  A  uniform  correction  has  been  determined  for  the  STD 
records,  and  is  described  in  paragraph  2  of  the  text,  but  has  not  been  applied  to  any  of  the  data  as  it  is  only  a  small  constant 
value.  All  temperatures  should  be  corrected  by  0.07°C,  all  salinities  by  -0.09°/ •». 
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from  2  to  20  km.  Vertical  scales  of  less  than  2  m  were  not  clearly  resolved  by  the  salinity,  temperature, 
and  depth  measurements  (STD)  actually  obtained.  However,  the  trace  pen  was  held  stationary  for  1- 
min  intervals,  and  the  heaving  of  the  ship  effectively  caused  the  tethered  sensor  to  sample  the  vertical 
temperature  structure  with  greater  resolution.  Figure  10  of  their  article  shows  this  sampling  at  a  275-m 
depth  (point  Pof  our  Fig.  9a,  their  Fig.  3).  From  this  figure  it  is  deduced  that: 

•  The  lowest  mean  temperature  was  17.26°C  at  minute  9  of  the  record. 

•  The  highest  mean  temperature  was  17.97°C  at  minute  39  . 

•  The  total  range  of  vertical  displacement  of  the  isotherm  due  to  internal  waves  was  5  meters. 

•  The  homogeneous  lamina  at  P  was  perhaps  3  m  thick. 

•  The  mean  temperature  gradient  at  station  13  between  200  and  300  m  was  4  x  10“ 2  °C/m. 

•  Inside  the  lamina  Pthe  gradient  was  less  than  3  x  10“ 3  °C/m. 

•  Directly  above  and  below  the  lamina  the  temperature  gradient  was  >  0.35°C/m  (3.5  x  10-3 
°C/cm). 

•  The  horizontal  scale  of  the  lamina  exceeded  2  km. 

Stommel  and  Fedorov  note  that  other  authors  (for  instance,  Liebermann  [12])  have  measured 
much  smaller  horizontal  scales.  They  conjecture  that  the  observed  lamina  of  their  own  measurements 
may  be  corrugated  so  that  a  horizontally  moving  probe  would  move  in  and  out  of  the  layer,  thus 
accounting  for  smaller  scales.  Liebermann  [12]  measured  temperature  fluctuations  in  the  ocean  along  a 
horizontal  line  by  means  of  a  probe  mounted  on  a  submarine.  His  work  is  reviewed  briefly  in  Appendix 
G.  Here  we  write  again  his  main  analytical  result,  that  the  mean  value  of  the  fluctuation  in  temperature 
gradient  is  given  by  the  formula 

<  [A  T{d)]2>  -  2<  T2>  [l  -  +  —  e~2  (units:  ( X:)2), 

in  which  the  unit  of  range  d  is  cm,  and  <  T2>  is  the  spatial  variance  of  the  fluctuating  temperature. 
The  horizontal  correlation  length  lc9  was  estimated  by  Liebermann  to  be  60  cm. 

It  is  somewhat  difficult  to  gather  from  this  work  what  the  value  of  <  T2>  should  be.  On  p.  565 
(of  Ref.  12)  there  appears  the  sentence,  "The  average  temperature  deviation  of  this  record  is  0.04°C.../' 
The  record  in  question  was  some  30  m  long,  taken  at  a  depth  of  —50  m  in  the  continental  waters  off 
the  U.S.  Actually,  a  variance  could  be  obtained  from  the  record  by  subtracting  out  the  mean,  selecting  a 
number  of  data  points,  squaring  the  temperature  fluctuation  at  these  data  points,  and  adding  all  the 
squares.  The  sum  would  be  the  variance.  The  square  root  of  the  sum  would  be  the  standard  deviation. 
This  has  not  been  done.  Since  we  desire  to  make  Liebermann's  data  compatible  with  GM72  [3],  we 
will  proceed  in  the  next  section  on  the  basis  of  this  latter  interpretation. 


The  Tow  Spectrum  of  Liebermann  as  Interpreted  by  GM72  (31 

Liebermann  [12]  measured  an  autocorrelation  of  temperature  fluctuation  of  form  —  e'p/60,  in 
which  p  is  the  log  distance  in  cm  and  60  cm  is  the  correlation  distance.  For  convenience  we  write  this 
as 

R  (p)  «  exp (— ti1  p)- 
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The  power  spectrum  of  these  temperature  fluctuations  can  be  obtained  by  use  of  the  Wiener-Kinchine 
theorem: 


FR(a)  cc  4  exp (~lcep)  cosap-^- 

a  1 

*  Uc~e2  +  a2) 


(for  derivation  see  p.  312  of  [9]).  Here  a  takes  on  positive  values  only.  The  constant  of  proportional¬ 
ity  is  the  mean-square  fluctuation  (or  variance)  of  the  random  temperature  field  for  all  spatial  and  tem¬ 
poral  microstructure  scales.  Thus  the  a -spectrum,  or  towed  spectrum  (TS),  of  Liebermann  as  inter¬ 
preted  by  GM72  [3]  has  the  form 


Ft  =  <h92>  — 

7 r 


IcV 


r2 


4*  or 


Now  at  fixed  depth  y  the  temperature  power  spectrum  is  related  to  the  vertical-displacement  power 
spectrum  F± : 


Fr  = 


d0_ 

dy 


h 


(see  p.  237  of  GM72  [3]).  Actually,  to  obtain  a  universal  representation,  GM72  plots  hF^{a\)  vsaj,  in 
which  the  following  definitions  hold: 


h  —  dimensional  Vaisala  frequency  (units:  cycle/h), 

^  (a  | )  =  dimensional  power  spectrum  for  vertical  displacement  (units:  m3/cycle),  and 
a  |  =  horizontal  component  of  wavenumber  (units:  cycle/m). 


To  check  this  representation  by  GM72  of  Liebermann's  work,  we  choose  a  depth  of  30  m,  such 
that  d0/dy  =  0.2  "G/m.  The  towed  spectrum  of  Liebermann’s  data  extends  roughly  from  0.1  to  10 
cycle/m.  GM72  interprets  Liebermann’s  measured  correlation  length  ( lc9  =  60  cm)  as«o  1  =  0(1  m). 
Inside  the  selected  wavenumber  band,  let  us  choose  aj  =  1  cpm  as  a  check  point,  and  let  us  take  h  = 
6  cycle/m  at  the  selected  depth.  Thus,  for  <S02>  —  (0.02)2, 


n  Fc(a i)  — 


(0.02)2 

(0.2)2 


2 
7 T 


1 

12+  l2 


x  6 


=  —  X  10-J  =  2  X  10-2  m3h_1. 

7 r 

This  checks  with  the  value  shown  in  Fig.  2  (Fig.  4  of  GM72).  If  we  were  to  construct  a  power  spec¬ 
trum  of  microtemperature  fluctuations  using  Liebermann’s  work,  we  could  use  this  result  to  find  one 
point  on  the  spectrum: 

Ft  =  (0.2)2  x  —  x  l<r2  x  =  —  x  l<r4 

7T  6  7T 

These  formulas  will  be  used  in  the  next  section. 


Calculation  of  Effect  of  Temperature  Microstructure  on  the  Phase  of  a  Fiber-Optic  Hydrophone 

As  noted  earlier,  this  report  uses  Liebermann’s  data  [12]  as  interpreted  by  GM72  [3]  to  find  the 
effect  of  temperature  microstructure  on  the  phase  stability  of  a  fiber-optic  hydrophone.  From  the  previ¬ 
ous  section. 
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Fr(a)  -  <802>  — 


/- 1 
lcU 


IcY  +  «2 


(units:  ( t)2  *  m). 


We  want  to  find  the  band-limited  variance  of  temperature  fluctuations: 


<802(&a)>  -  f  <S02>  — 

Aar  7T 


/-I 

Irti 


iiY  + 


c/a  (units:  (TT)2). 


The  integral  is  elementary: 

<802(Aa)>  =  <802>  |arctan(amax/t^)  -  arctan(amin/^)J. 

Now  choose  a  fiber  L  units  long  and  assume  that  the  fiber  spatially  detects  the  temperature  fluctuations 
along  its  entire  length.  Then  the  bandwidth  of  wavenumbers  extends  from  °°  to  ML.  Thus, 


<802(Aa)>  =  <802>  — 


For  example,  let  L  —  40  cm  and  Ic9  =  60  cm;  then 


<802(Aa)>  -  (0.02)2  — 


t  #  t# 

arctan  oo  —  arctan— 


60 


7  "  arCt3n  40 


Similarly,  for  L  —  1  cm, 


-  1.50  x  10“4  (units:  (X’)2). 


<802(Aa)  >  -  4.24  x  10’6  (units:  (V)2). 


The  effect  of  these  temperature  fluctuations  on  the  microstructure  scale  is  found  from  the  for¬ 
mula  previously  derived  (in  the  section  "Calculation  of  the  Mean-Square  Fluctuation  in  Temperature 
over  a  Towed  Fiber"): 

V<802>  =  1.400  x  102  V<802(AoT)>. 

Thus,  the  mean-square  phase  fluctuation  caused  by  temperature  microstructure  is  calculated  to  be 
V<8</>2>  =  1.400  x  102  x  Vl.5  x  10“4  =  1.7  rad,  for  L  =  40  cm. 


and 


V <8<£2>  =  1.400  x  102  x  V4.2Tx  10  6  =  0.29  rad,  for  L  =  1  cm. 

The  magnitudes  of  these  phase  disturbances  are  contingent  on  the  selection  by  GM72  [3]  of  dT  = 
0.02  °C  and  8T/8y  =  0.2°C/m.  A  nearly  similar  result  is  obtained  if  we  use  Liebermann’s  measure¬ 
ments  directly  1 12].  As  explained  in  Appendix  G,  Liebermann  measured  an  average  fluctuation  of 
about  0.04°C  over  a  correlation  distance  of  about  60  cm.  By  extrapolation  the  rms  average  fluctuation 
over  a  distance  of  1  cm  is  about  3.4  x  10_3oC,  and  over  a  distance  of  40  cm  it  is  3.9  x  10_2oC.  Given 
the  uncertainties  of  the  precision  of  measurement,  these  values  are  near  enough  to  make  the  interpreta¬ 
tion  of  Liebermann's  work  by  GM72  a  reasonable  model. 

CONCLUSIONS 

The  temperature  field  in  the  ocean  is  subject  to  spatial  and  temporal  fluctuations  of  two  distinctly 
different  scales.  The  larger  scales,  of  the  order  of  hundreds  of  meters  in  space  and  fractions  to  several 
hours  in  time,  are  caused  by  internal  waves.  The  smaller  scales  (microscales),  of  the  order  of  millime¬ 
ters  in  space  and  seconds  in  time,  are  caused  by  velocity  turbulence  and  thermal  diffusion.  We  have 
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calculated  the  effect  of  both  these  types  of  temperature  fluctuations  on  the  phase  stability  of  an  optical- 
fiber  hydrophone  when  the  summation  is  nearly  incoherent.  The  results  of  the  analysis  appear  in  the 
form  of  space  (subscript  L)  and  time  (subscript  7)  phase  fluctuations: 

ceLy/2LTcserms(L). 

and 

(8</>  7~)rms  =  C»T  yfeTTc  $Qrms(T), 

in  which  (8 9L)Tms  —  V<802(Aa)>  and  (80f)rms  *=  V <802(Aa>)> ,  where  <802>  is  the  mean-squared 
fluctuation  of  temperature  over  either  a  band  Aa  of  wavenumbers  or  a  band  Aco  of  frequencies.  To 
obtain  the  latter  we  have  used  the  work  of  Garrett  and  Munk  and  of  Liebermann.  Table  3  shows  80rms 
calculated  for  specific  choices  of  spectrum,  fiber  length,  and  depth  of  ocean,  based  upon  the  assumption 
that  \2LLC\  =  0(1)  m2  and  |27Tf  |  -  0(1)  s2. 


Table  3— Calculated  RMS  Temperature  Fluctuations  in  the  Ocean 
and  Corresponding  RMS  Phase  Fluctuations 


Method 

of 

Analysis 

Literature  Source 
of  Spectrum 
or  Correlation 

Length 

of 

Fiber 

(m) 

Depth  of 
Center  of 
Fiber  in 

Ocean 

(m) 

RMS  Vertical 
Displacement 
of  Ocean  in 
Spectral  Window 

(m) 

RMS  Temperature 
Fluctuation  in 
Spectral  Window 

(°C) 

RMS  Phase 
Fluctuation  for 
\2LLC\  -  Od)m2 
or 

|27Trl  -  0(l)s2 
(rad) 

Power  spectrum  of 

internal  waves 

Horizontal 

GM72  [3] 

1000 

100 

1.15 

0.029 

4.04 

(Towed) 

Vertical 

GM72  (31 

100 

100 

1.1 

0.027 

3.83 

(Dropped) 

Temporal 

GM72  131 

100 

100 

5.7* 

0.14* 

19.6* 

(Moored) 

1.13+ 

0.028+ 

3.92+ 

0.46* 

0.0114* 

1.6* 

Spatial  correlation  of 

Liebermann  (121 

0.4 

30  to  60 

_ 

0.012 

1.7 

temperature  microstructure 

(Correlation  length  «  60  cm) 

0.01 

30  to  60 

- 

0.0021 

0.29 

‘"mil/"/  ~  1 
+"tnir/"/  ”  25 


Aw  band  “* 


"max  *  3  cPh 
"mm"  0049  cPh 


*"mir/",  “  150 


The  significant  temporal  contribution  to  the  power  spectrum  of  internal  waves  is  contained  in  the 
band  between  the  inertial  frequency  (0.04  cph  in  our  latitude)  and  the  Vaisala  frequency  (3  cph  and 
upward  depending  on  depth).  The  frequencies  of  interest  for  a  fiber-optic  hydrophone  are  currently 
above  10  Hz;  therefore,  it  is  concluded  that  the  temperature  fluctuations  in  time  associated  with  inter¬ 
nal  waves  will  not  directly  contribute  to  fiber-optic  hydrophone  noise.  Additional  noise  will  not  be 
present  even  if  a  fiber-optic  sensor  is  towed  at  normal  speeds  through  the  temperature  of  internal  waves 
on  a  level  track. 

Lack  of  data  has  prevented  the  determination  of  temporal  dependence  of  temperature  microstruc¬ 
ture;  however,  spatial  calculations  indicate  that  a  reduction  of  unwanted  phase  shifts  due  to  temperature 
microstructure  can  be  achieved  by  the  use  of  an  extended  sensor  (see  Table  3,  "Spatial  correlation  of 
temperature  microstructure").  When  additional  data  concerning  the  relative  temperature-induced  phase 
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shifts  of  piezoelectric  and  fiber-optic  sensors  are  available,  analysis  of  the  temporal  nature  of  tempera¬ 
ture  microstructure  may  be  valuable. 

The  rms  phase  fluctuations  due  to  temperature  effects  of  both  internal  wave  and  micro¬ 
structure  indicate  that  homodyne  fading  could  be  significant  as  the  operating  point  is  shifted  away  from 
the  optimum  phase  setting.  The  magnitude  of  these  fluctuations  will  be  important  if  homodyne  detec¬ 
tion  with  feedback  for  optical  alignment  is  used.  The  frequency  response  of  a  feedback  system  would 
depend  on  the  temporal  dependence  of  internal  waves  and  microstructure. 

It  is  important  to  note  that  the  results  listed  in  Table  3  are  strongly  depth  dependent  and  that  they 
are  narrowly  contingent  on  the  use  of  band-limited  spectra.  Furthermore,  we  have  assumed  that  the 
temperature  fluctuations  are  instantly  sensed  by  the  optical  fiber,  and  we  have  thereby  neglected  the 
finite  time  effects  of  thermal  conduction.  Formulas  which  include  conduction  have  been,  however, 
explicitly  derived  in  this  report.  It  is  natural  to  expect  that  thermally  insulated  fibers  will  exhibit  greatly 
reduced  sensitivity  to  ocean  temperature  fluctuations. 

In  a  separate  analysis  we  have  derived  formulas  for  the  S/ N  ratio  of  fiber-optic  hydrophones  in 
which  the  noise  V  includes  laser  jitter,  environmental  noise,  and  the  customary  shot  noise. 
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Appendix  A 


CORRELATION  FACTORS 

From  the  "Procedure  of  Analysis"  section  we  can  find  that 


<S<f>l>  =  2C2  L  J0 1 


Ru  (f)^. 


where  the  units  of  C  are  rad/°C  •  m,  of  L  are  m,  of  R ^  are  (°C)2,  and  of  f  are  m. 
Now  assume  that  is  Gaussian, 


Rm  =  <802>  e 


2^ 


then 


Now  [A  1  ] 


<ul>  =  2 C2L  <802>  j/o"  e~ilnL<dt  -  // "f;  ^ 

Jo 


di 


•l  -f}/2L,2  „  yfrr  ^ 

e  di - ; - r  <J> 


V2Lf 


V2Z., 


<J>lx)  =  -t=-  f  ?  '2 

V7T 


and 


xJo  ^ —  (i - e  >• 


Thus, 


<h<t>l>  =  2C2L  <802> 


When  Lt  «  L,  the  correlation  factor  is  Lc  V?  <t>  (00)  —  Lc  V? 


The  function  1  -  acts  as  a  "window"  to  the  correlation  function.  Assume  that 


periodic,  say  <802>  cos  aL  Then 


'8#^)  is 


s  nanism 


_ 


If  —  «  1,  then  the  correlation  length  is  very  large  and  the  correlation  factor  K  is  L/2.  If 

kL/2  »  1,  then  there  are  numerous  oscillations,  the  correlation  length  is  very  small,  and  the  correla¬ 
tion  factor  K  — >  0. 
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Appendix  B 

MEASURED  VS  THEORETICAL  AUTOCORRELATIONS 


A  review  of  an  analysis  by  Corcos  [Bl]  is  presented  here. 


Let  Nix,  t )  be  the  instantaneous  disturbance  of  the  medium  at  point  X,  and  let  < N(x)>  be  its 
average  value  over  a  long  duration.  The  fluctuation  of  N  is  defined  as  nix,  t)  =  Nix,  /)  -  <  Nix)>. 
It  is  considered  a  random  variable  in  space  and  time  and  possesses  an  autocorrelation  function 

R  (£,  t)  =  <nix,  t)  nix  +  £,  t  4-  r)> 

and  a  power  spectrum  (or  energy  density) 

oo 

Fia,  oj)  ~  f  f  f  R  (£ ,  r)  exp[+  /(cor  -  a  •  £)]  . 

•'r0r  (2irr 

The  partial  Fourier  transform,  or  cross-spectral  density,  is  defined  as, 

F(£,a>)  =  f  *(£,  r)  e+,u>T-~~. 

Jo°  2n 

Now  consider  the  measurement  instrument  to  be  a  filter  with  weighting  function  (or  window)  his,  17) 
(units:  m~3s“!).  The  measured  disturbance  is  then  the  convolution, 

00 

nix,  t)  =  f  f  f  f  n(x>  f)  Ms  -  x,  7)  /)  ds}  dr). 


The  measured  autocorrelation  over  volume  Fand  time  Tis 

Rm  =  §  §  •  •  •  J*Ms-x,  17—  /)  Mx,  /)  Mx  +  £,  t  +t)  //[q  -  (x  +  £),f  -  (/  4-  rl 

Let  s  —  x  «  a,  q  -  (x  +  £)  =  b,  17  -  /  =  c,  and  f  —  (t  +  r)  =  d,  and  consider  spatial  coordinates  first. 
Then  for  stationary  homogeneous  statistics  of  n  the  integrand  has  the  form 

//  (a)  n  is  -  a)  n  (q  -  b)  h  (b) . 


The  total  spatial  lag  is 

s-a-q  +  b=S“q  +  b-  a  =  {+  b-a. 

Hence  an  integration  over  d}x  yields  the  autocorrelation  in  space 

Rii  +  b  -  a). 

A  similar  consideration  holds  for  the  temporal  coordinate  so  that  an  integration  over  dt  yields 

Rir  +  d- c). 

Thus  the  measured  autocorrelation  function  is 


00 

RJt.  T)  =  JJ  -  }h(a.c)R{l+b-»,T  +  d-c)h(b,  d)  d'td'bdcdd. 


47 


S.  HANISH 


This  form  allows  us  to  derive  the  measured  cross-spectral  density  at  zero  lag  (£  -  0),  which  is  the  ordi¬ 
nary  frequency  spectrum 4>m(<*>): 

<Dm  M  =  F„,  (0,  (o  >  =  Rm  (0,  t  )  e'“'  jL 

oo 

—  J  e,liiT~-  J  J  *  c)R  (b  —  a,  r  +  d  ~  c)  h  (b,  d)  dP a  c/3b  dc  dd. 

—  oo 

An  oft-encountered  case  is  the  one  in  which  the  transducer  instantaneously  senses  the  effect  of  the 
environmental  disturbance.  In  this  case 

hi a,  c)  Mb,  d)  =  hi a)  Mb)8(r  -  d). 


When  this  applies,  the  following  formulas  hold: 


=  S-oo  e *  J*  (b—  a,  T)h(b)dySkdyb 

— oo 
oo 

-//  ■  J  //(a)  Mb  -  a,  ct>)  //(b)  c/3a  c/3b. 


Now  let 


b  —  a  =  e, 


then 


4>m(o>)  -  ///««.  co)  //(a)  /? (a  +  e)  d} a c/3e. 


From  this  one  can  identify  a  transducer  response  function  H(e), 

//(e)  =  f  /i(a)/i(a  +  e)rf3a. 

*/  oo 

This  is  sometimes  called  the  autocorrelation  function  of  the  filter  window.  The  measured  frequency 
spectral  density  is  then 

oo 

(to )  —  us  Fie,  to)  //(e)  d3e. 

*—  OO 

Now  Fie,  o) )  is  related  to  the  power  spectrum  of  the  total  (temporal-spatial)  disturbance,  I  Mia,  oj)|2: 

Me,  co)  =  f  |  Alia,  oj)  |2  e~a  e  d}a , 

where 

A  more  general  transducer  response  function  occurs  when  the  temporal  response  is  finite  in  time.  Let 

d  —  c  —  g. 


so  that 


Hie,  g)  =  ///«..  c)  h  (a  +  e,  c  +  g)  dy a  dc. 

—  oo 


The  measured  temporal  spectral  density  is  then 


Oo 

d>m  (o»)  =  us  Fie,  g,  co)  hie,  g)d}edg, 

—  oo 
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where 


F(e,  g,u>)  =  J  f  I M(a2.  w)|2  e~‘a  e  e+i““  dvd'*. 

—  oo 


Window  effects  may  be  approached  in  another  way.  Consider  spatial  effects  alone.  Let  w(x)  be  a 
spatial  window  (units:  none).  Then 


W(k)  =  r  e~iK  x 

— oo 


(units:  m3). 


If  F(y)  is  a  Fourier  power  spectral  density  of  an  environmental  disturbance,  it  can  be  "seen"  through 
the  instrumental  window  w(x): 

F(a)  =  J  w  W(a  —  y)  F( y)  d*y. 

Example:  Let  the  window  be  a  triangle  centered  on  x  of  base  2 L\ 
w(A-)  =  i  -  lx |  <  L, 


Then 


w(x)  =  0 

-ii‘ 

7T  JO 

^  JL 

2tt 


1*1 

L 


1  - 


L 

sin(*L/2) 


otherwise. 

„f «V  dx_ 

27 T 

COSK.V  dx 
2 


K  L 
2 


(units:  m). 


If  F(y)  is  the  power  spectral  density  of  a  disturbance,  then  the  convolution  gives  the  " windowed"  spec¬ 
trum: 

12 


L  r~  lsini (y_a) 


(y  —  a) 


F(y)  dy. 


define  a  band  of  wavenumbers  Ac*  by  the  requirement 

AaJ_, 


2tt 


The  integral  is  equal  to  7r.  Therefore  the  band  of  wavenumbers  is 


r 

.  L 

Siny  a 

2 

i  L 

J  —  oo 

L 

~a  | 

C 

Ac*  ~  y 


(units:  m  ’). 


in  which  L  is  half  of  the  base  of  the  triangle  representing  the  window.  The  effective  band  (at  the  3-dB 
power  points)  is  half  of  this.  Thus, 

L 


A  a 


e  flee  live 
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Appendix  C 

FORMS  OF  FOURIER  INTEGRALS  USED  IN  THIS  REPORT 


Let  lit)  be  a  random  current.  Then  it  can  be  represented  as  a  Fourier  expansion, 

/(/)  =  r  e-ia“l(<o)dto, 

%f  —  oo 

in  which  /(o>)  is  the  amplitude  spectrum,  with  units  A  s.  By  inverse  expansion  one  recovers  /(/): 

/(/)  =  f  °°  e~iul  f°°  i(t')  e+iu'  4L  d(u 

This  is  formally  proved  by  identifying  the  delta  function, 

8(/  -  r')  -  ~  r 

2n  J  ” 


In  this  report  we  associate  the  2n  in  the  expansion  of  /(/)  arbitrarily  with  the  time  coordinate, 
i.e., 

/(«)  =  /”  tu) 

•/-oo  Z7T 

Particularly,  in  the  application  of  Rayleigh’s  theorem,  we  use  the  convention, 

J-f  \i(t)\2dt—  f  |/(co)|2  da). 

2.7T  **  ~°° 

The  units  of  the  power  spectrum  |/(o>)|2  are  A2  s2.  If  the  time  duration  of  /(/)  is  finite,  say  T s,  then 


f  \i(t)\2  dt 

J  —  oo 


T  </2>, 


in  which  </2>  is  the  mean-square  amplitude  over  T.  Thus  for  finite  signals  (in  time)  Rayleigh’s 
theorem  has  the  form 

I  K(*))\2 


<il> 


do), 


in  which 


Thus, 


/7Y2  dt 

_  /(/)  e,w(  — — . 


-77  2 


27T 


2tt  [/  (a>)  1 2 

T 

is  the  power  spectral  density  (in  units  of  A2  •  s)  of  the  random  function  /(/)  in  the  interval 
-772  <  t  <  77 2,  and  /(w)  is  the  Fourier  transform  of  /(f). 

These  same  conventions  are  applicable  to  spatial  as  well  as  temporal  Fourier  transforms.  If 
p( r,  /)  is  a  random  function  of  space  and  time  by  temperature  (units:  °C),  then  the  mean-square  p  over 
all  space  V  and  time  T  is 

<p2>  -=  (2rr)4  JJ  .If  <K^IL  dKdo)  (units^TT)2), 


S.  NANISM 


where 


/>(K,w)  =  f  f  p(r,  t)e  'K  -~~r  (units: 'C ■  m3 •  s). 

J  J  2n  (2 ttY 


Thus 


(2tt)4Ip (K>  u))\2 
VT 


(unitsiCC)2  *  m3  *s) 


is  the  power  spectral  density  of  the  random  function  p( r,  /).  Similarly,  if  the  spatial  part  of  the  spec 
trum  is  one-dimensional,  then  the  power  spectral  density  is 

(2tt)2  (units:(‘C)2*  m  s), 

where 

p(a,(o)=  f  p(x,  t)  etwl  ~  (units:  <C*m*s) 

J-oo  (2tt)  2 

and 

<p2>  »  f  (27r)2  dado)  (units:  (XI)2). 


Appendix  D 


POWER  SPECTRAL  DENSITY  FACTORS  FOR  sinx/x  CORRELATION 
OF  THE  RANDOM  FIELD 

Let  ij(r )  be  the  random  current  caused  by  laser  jitter.  Assume  that  the  autocorrelation  function 
of  ij(t)  has  a  sinx/x  dependency.  Then 

sin7r  (QJ<a)ct) 


Y(t)  =  <  ij  > 


7t(0Jojct) 


(units:  A2). 


This  form  is  selected  because  it  closely  resembles  the  time-delay  correlograms  of  ambient  ocean  noise 
obtained  by  Urick  [Dl].  In  addition,  when  q)ct  =  1, 

sin7r  (0.7) 


tt(0.7) 

showing  that  the  autocorrelation  cuts  off  as  required.  Thus  tc  becomes  the  correlation  time.  Again, 
according  to  the  Wiener-Kinchine  theorem. 


j_r 

In  J-° 


e~tar  Y (r)dT  =  y  |S5y(w)|2 


(units:  A2  s). 


Now 


where 


and 


l  r°° 

-i-j  e~iQtT\{r)dr 

ITT  J~°° 


B 


0.7 (oc  x  2t r  7rx 


smrxe-l2nXSdx^ 


B 


0.7ct>f  x  27t 


n  (s). 


x  =  0.7cocr 
B  =  </7>, 


5  — 


to 


and 


0.7  x  2ntoc 
n(s)  is  the  rectangular  function 

U(s)  -  1  when  s  = 


r  _  -  <  4*,  or  co  <  0Jntoc, 

0.7  x  27tojc  2  c 

l 


n(s)  -  0  when  s  >  y,  oro>  >  0.l7Ttoc. 
Thus,  the  laser  jitter  power  snectral  density  factor  is  modeled  as 


2n 


1 8  Sj  (co )  | 2 


<02>n 


to 


0.7  x2ntor 


T  0.1toc  x  2t r 

The  autocorrelation  and  nower  spectral  density  factor  are  shown  in  Fig.  Dl. 

A  similar  form  can  be  obtained  for  one-dimensional  spatial  autocorrelations,  viz.. 


Y  (/)  — 


:/f2f>  sinTr  0.7-/“ 
Ic 


(units:  A2). 
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Fig.  D1  —  (a)  Autocorrelation;  (b)  Power  spectral  density  factor 


The  power  spectral  density  factor  over  both  space  and  time  is,  therefore. 


(2„); 


aic 


1 


0.7  x  2tt 

[ 

n 


0.7wr  x  2rr  J0.7x27rc*jc 
in  which  L  is  the  spatial  interval  of  integration. 


(units:  A2  ■  s), 


Summary 

The  noise  contributions  (units:  A2/ Hz)  from  fluctuations  in  a  laser  beam  and  from  environmental 
effects  have  been  modeled  on  the  basis  of  sinx/x  autocorrelations.  The  resulting  formulas  for  the 
noise  are  written  in  terms  of  rectangle  functions  II  (s)  *=  US  ^  1/2,  and  IlCs)  =  0,  S  >  1/2 


Laser  beam  jitter:  The  corresponding  random  photomultiplier  current  is  /)(/),  whose  power  spec¬ 
tral  density  factor  is  modeled  as 

<i}> 


0.7*2™, 


n 


0.7  x  2tT(Ocj 


(units:  A2-s), 


where  u>cJ  is  the  correlation  frequency. 


Space -  and  time-varying  environmental  temperature :  The  corresponding  random  photomultiplier 
current  is  /w( r,  r),  whose  power  spectral  density  factor  is  modeled  as 


(27t)2 

LT 


|6SW  (a,o>) 


(2tt)2x(0.7)V, 


n 


ai 


\ c9 


n 


0.7  x  2tt  I  0.7  x  2ir  c*ce 


where  is  the  correlation  frequency,  and  lc8  is  the  correlation  length. 


3-space  and  time-varying  environmental  mechanical  stresses:  The  corresponding  random  photomulti¬ 
plier  current  is  /w(r,  t),  whose  power  spectral  density  factor  in  an  isotropic  environment  is  modeled  as 
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2tt 

Vp 


IS*Sa/(/4, 


<o)\2  = 


<ti>  (rA 


M 


(2tt)4x(0.7)4o>cM 


n 


M 


0.7  x  2n 


II 


a> 


0.7  x  27 to)cM 


(units:  A2  m3  s) 


where  (ocM  is  the  correlation  frequency,  lcM  is  the  correlation  length,  and  p.  -  |jt  |. 
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Appendix  E 

EQUIVALENT  CIRCUITS  OF  PHOTODETECTOR  AND  PREAMPLIFIER 

A  generalized  pholodetector  can  be  schematically  represented  as  an  input-signal  current  source 
'signal  in  parallel  with  an  admittance  Yf  and  a  noise  current  source  /£.  The  internal  current  gain  is 
<  G  > .  The  circuit  is  shown  in  Fig.  El . 


The  noise  current  source  is  typically  shot  noise.  Using  symbols  defined  in  the  text,  this  noise  is 
<iS>  =  2elUph)toh  +  (/,*)incoh  +  lBK  +  ld]  <G>2  F(G) 
in  units  of  A2/Hz.  It  is  to  be  noted  that  /sjgnal  =  Uph)m h  +  (^)incoh  and  that 

Y/  =  Rf'  +  j<D  Q. 

A  generalized  preamplifier  can  be  schematically  represented  as  a  noiseless  generator  with  gain 
A  (w)  and  input  admittance  YA  =  GA  +  j(*>CA .  It  has  two  noise  sources  [El),  a  voltage  noise  source  vA 
in  series  with  the  input  terminals  and  an  uncorrelated  current  noise  source  iA  in  parallel  with  the  input 
terminals.  Coupled  to  iA  is  a  second  current  source  vA  Yc,  which  is  correlated  to  v A  through  a  cross 
admittance  Yc  ==  ReFc.  +  j\m  Yc.  The  equivalent  circuit  is  shown  in  Fig.  E2. 


vA 
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Appendix  F 

NUMERICAL  CALCULATIONS  OF  MEAN-SQUARE  PHASE  FLUCTUATIONS 
CAUSED  BY  TEMPERATURE  FLUCTUATIONS  IN  THE  OCEAN 

The  formula  derived  in  the  text  is 

V<80/>  =  1.40  x  102  V<802(Aa,  Aw)>. 

As  noted  in  the  section  "Calculation  of  the  Mean-Square  Fluctuation  in  Temperature  over  a  Towed 
Fiber,"  the  mean-square  temperature  fluctuation  for  a  particularly  specified  fiber  1000  m  long  at  a  depth 
of  100  m  is 

<802  (Aa)>  -  8.31  x  10~4  (°C)2, 
and 

V<8<*)2>  =  1.40  x  102  x  V8.31  x  10~4  =  4.04  rad. 

Calculations  for  L  —  100  m,  10  m,  and  1  m  are  shown  below. 

L  **  100  m,  a  j  =  10“ 2  cpm: 

1  min 

/  =  (1.22  x  10“4)2  x  — =  1.4884  x  lO-6, 

1(T2 

f  Fc  (a,)  dax  =  8.925  x  104  x  1.4884  x  10^ 

=  1.33  x  10"'  m2, 

<S02>  =  (0.025)2  x  1.33  x  10“ 1  =  8.31  x  10~5  (°C)2. 

and 

A<£  =  1.40  x  102  x  V8.31  x  10“5  =  1.28  rad. 


L  «  10  m,  a.  =  1(T 1  cpm: 

‘min 

/  =  (1.22  x  10“4)2  — p-  -  1.4884  x  10"7, 

f  Ft  (a,)  da,  -  8.925  x  104  x  1.4884  x  10~7 
=  1.33  X  ir2  m2, 

<S02>  =  (0.025)2  x  1.33  x  10~2  -  8.31  x  10“6, 
and 

A<£  -  1.40  x  102  x  VOI  x  10~6  -  0.404  rad. 


L  *  1  m,  a  |  =1  cpm: 

1  min  r 

/  =  (1.22  x  10“4)2  x|-  1.4884  x  lO’8, 
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Appendix  G 

STRUCTURE  FUNCTIONS  OF  TEMPERATURE  AND  LIEBERMANN’S  MODEL 


Let  C(x,  t)  be  the  speed  of  sound  at  point  x  in  the  ocean.  Assume  that  the  ocean  is  "frozen" 
temporarily  and  plot  C(x)  vs  x.  It  is  found  that  C(x)  is  a  random  function  with  mean  C0(x)  and  a 
fluctuation.  Suppose  that  the  mean  value  is  subtracted  out,  and  designate  the  remaining  fluctuating  part 
by  8C(x).  Now  at  point  x  one  can  find  8C2(x).  The  ensemble  average  of  8  C7(x)  for  all  x  is  its  vari¬ 
ance.  The  square  root  of  the  variance  is  the  standard  deviation  of  C(x)  at  x,  i.e.. 


f <  C2(x)> . 


Between  two  points,  x  and  x  +  p,  there  will  be  a  finite  difference  in  C(x).  By  selecting  a  large  number 
of  two-point  pairs,  one  can  construct  the  ensemble  average,  or  structure  function,  Dc  (p), 

A-(p)=<l8C(x  +p)  — 8C(x)P>  =  <8C2(x  T  p)>  +  <SC2(x)>  -  2<8C(x  +p8C(x)>. 

Assume  that  8C(x)  is  statistically  stationary  in  space,  then 

<8C2(x)>  =  <  8  C2  (x  4-  p )  >  =  <8C2(0)>. 

Therefore, 

Dc(p)  =  2[RC(0)  -  Rc(p)]t 


in  which 


Rc(p)  =  <8C(x  4-  p)8C(x)> , 

i.e.,  Rc(p)  is  the  autocorrelation  of  the  zero  mean  fluctuation.  Dividing  through  by  Rc{ 0),  one  has 

<[8C(x  +p)-SCU)l2>  =  J.  _  Ryip) 

<8C2(0)>  [  RAO)  ' 

This  is  the  average  difference  squared  between  two  points  of  the  frozen  fluctuation  pattern  normalized 
to  the  variance  at  a  point.  Since 

C(x4-p)-C(x)==  C0(x  4*  p)  +  8C(x  +  p)  —  C0(x)  —  8C(x), 

we  can  say  that,  for  a  short  distance,  C0(x+p)-  Co(x)  =  0.  Hence,  we  can  say  that 
8C(x  4-  p)  —  8C(x)  measures  the  spatial  gradient  of  the  speed  of  sound. 


The  speed  of  sound  is  a  function  of  temperature,  salinity,  pressure,  etc.  In  the  first  approxima 
tion,  for  temperature  J(°C),  assume 

C(x)  —  C0  +  421 3Tx)  4-  . . .  (units:  cm/s). 

Thus,  if  8  Tlx)  is  the  zero-mean  fluctuation  of  J,  then  AC(x)  «  421  8  T(x)  and 

C (x  +  p )  -  C(x)  -  421[8r(x  +p)  -  6T(x)l 

Since  we  are  normalizing  with  respect  to  the  variance  of  the  fluctuation  in  C(x)  at  a  point,  we  see  that 


<  (8T(x  4-  p)  -  8T(x)]2> 

<8  T2(0)> 


Rr^p) 
Rt( 0)  * 


Liebermann  [Gl)  measured  temperature  fluctuations  in  the  ocean  at  depths  of  30  to  60  m.  He 
found  that  the  mean  temperature  fluctuation  at  a  point  averaged  over  his  record  was  about  0.04°C  and 
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that  the  correlation  length  in  temperature  fluctuation  between  points  was  about  60  cm.  The  mean-square 
fluctuation  in  index  of  refraction  was  5  x  10“9.  Thus,  between  any  two  points  along  a  track  at  these 
depths  the  ensemble-averaged  difference  in  temperature,  i.e.,  difference  in  fluctuation  of  temperature 
(assuming  the  mean  is  subtracted  out),  is 

<[ST(x  +p)  -  8T(j c)]2>  =  <6T2  (p)>  =  2  x  0.042  [1  -  RT(p)]t 
in  which  R^ip)  is  the  normalized  autocorrelation  of  the  zero-mean  function  bT(x). 


Now  Liebermann  also 
R(p)  = 


found  that  the  normalized  autocorrelation  in  C(x)  can  be  modeled 
30  _ p  j  |  ^ 

^  e  60  ~  e  2  >  P  =  separation  distance  of  points. 


as 


Hence  the  horizontal  temperature  gradients  in  the  ocean  due  to  fluctuations  (with  mean  values  sub¬ 
tracted  out)  are  given  by 


V<A  T2(p)>  -  J2  x  0.04 


30e-&  +  _L/f 

29  29 


1/2 


(units:  X). 


Temperature  gradients  based  on  this  formula  are  given  in  Table  Gl. 


Table  Gl  —  Temperature  Gradients 


p 

1  -  Rip) 

V<AT2(p)> 

(cm) 

(°C) 

0.1 

4.0957  x  10- 5 

3.62  x 

10-4 

1.0 

3.53058  x  103 

3.36  x 

10- 3 

40 

4.6887  x  10  ' 

3.87  x 

10-2 

60 

6.194  x  10“' 

4.45  x 

10“2 

oo 

1 

5.66  x 

10-2 
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